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1, With rddpciot to ^^ionia liquid" 

The term "'ionic liquid" is the general terms (shown 
in References 1 to 5 as attached) . 

<ReferenGe 1> 

Masahiro Yoshizawa, Wataru Ogihara, and Hiroyuki Ohno, "'Novel 
Polymer Electrolytes Prepared by Copolymerization of Ionic 
Liquid Monomers"/ Polymers for Advanced Technologies, 13, p.589- 

594 (2002) 
<Reference 2> 

Masafumi Yoshio, Tomohiro .MuJcai, Kiyoshi Kanie, Masahiro 
Yoshizawa, Hiroyuki Ohno, and Takashi Kate, '''Liquid-Crystalline 
Assembiiea Containing Ionic Liquids: An Approach to Anisotropic 
Ionic Materials", Chemistry Letters 320-321 (2002) 
<Reference 3> 

Ann E. Visser, Matthew Reichert, Richard P. Swatloski, 
Heather D. Willauer, Jonathan G. Huddleston, and Robin D. Rogers, 
""Characterization of Hydrophilic and Hydrophobic Ionic Liquids: 
Alternatives to Volatile Organic ' Compounds for Liquid-Liquid 
Separations", ACS Symposium Series (American Chemical Society) , 
p. 289-308 (2002) 
<Reference 4> 

Masayoshi Watanabe, Tomoo Mizumura, ""Conductivity study on ionic 
liquid/polymer complexes". Solid state Ionics, p. 353-356 (1996) 
<Reference 5> 

Timothy I. Morrow and Edward J. Maginn, ""Molecular Dynamics 
study of the Ionic Liquid l-n-Butyl-3-methylimida20iium 
Hexafluorophosphate", Journal of Physical Chemistry B, 10 6,. 
p. 12807-12813 (2002) 

The term ""ionic liguid^^ refers to a salt which is in 
the liquid state at room temperature (may be called 
""ambient temperature molten salt"), and particularly to a 
salt having a melting point of 70**C or less, preferably 
30**C or less. Such an ionic liquid has no vapor pressure 
(non-volatility) , exhibits high thermal resistance and 



MUTUAL NON-DISCLOSURE AGREEMENT 



This Mutual Non-Disclosure Agreement ("Agreement") effective as of the last date of the signatures hereto (the 
"Effective Date") between RF Micro Devices, Inc. and its wholly owned subsidiaries ("RFMD"), a corporation having its 
principal place of business at 7628 Thomdike Road, Greensboro, North Carolina 27409, and Diemat Incorporated, a 
corporation having its principaLplace of business at 19 Central Street, Byfield, MA 01922 ('^PARTICIPANT') 
(individually a "Party" and collectively the "Parties"). 

WHEREAS, the Parties desire to enter into discussions regarding technology for mixing offerrite paste and its 
potential applications (the 'Turpose"); and 

WHEREAS, the Parties desire to disclose certain Confidential Information (as defined below) on a confidential 
basis to fiuther the Purpose. 

NOW THEREFORE, in consideration of the disclosure of Confidential Information by either Party, the Parties 
agree as follows: 

1. Definitions: 

a. The Party receiving Confid^tial Information firom the other Party shall be the "Receiving Party" 
and the Party disclosing Confidential Information to the other Party shall be the '"Disclosing Party." 

b. "Confidential Information" shall mean any information or data, whether owned by a Party or its 
AffiUates or obtained from a third party subject to a confidential arrangement, including, but not limited to, 
'Trade Secret Information" and "Business Information," ^sclosed by one Party to the other, in coimection with 
the Purpose, irrespective of the medium in Wiich such information or data is embedded, which infijrmation if in 
tangible or electronic form is marked as proprietary or confidential by the Disclosing Party, is 

designated by the Disclosing Party as being confidential at the time of such disclosure; oriwhich under the 
drcumstances surrounding tlie diisdosure gives the indicia that the information should be treated as confidmtial. 
Confidential Information shall include any copies or abstracts made thereof as well as any modules, samples, 
prototypes, or parts thereof embodying such Confidential Information. 

c. 'Trade Secret Information" includes, but is not limited to, product specifications, data, technical 
know-how, formulae, processes, designs, sketches, photographs, graphs, drawings, samples, prototypes, 
inventions and ideas, past, current, and plaimed research and development, current and planned manu&cturing 
methods and processes, computer software and programs (including object code and source code), draft or 
pending patent s^plications, and other similar technical materials, however documented. 

d. '"Business Information" includes, but is not limited to, information concerning the business and 
afl&irs of the Disclosing Party, which may include unpublished financial projections and budgets, projected sales, 
capital spending budgets and plans, persoimel information of the Disclosing Party, marketing plans, and other 
business information. 

e. "AflBliate" shall mean any company or other entity of which, and as long as, at least fifty percent 
(50%) of the outstanding shares or securities representing the right to vote for the election of directors or other 
managing authority are, now or hereafter, directly or indirectiy owned by a Party. 

2. Confidential Information of the Disclosing Party shall remain the property of the Disclosing Party. The 
Receiving Party agrees to protect the Confidential Information of the Disclosing Party against unauthorized disclosure 
with at a minimiun the same degree of care as it appUes to its own confidential information of like nature, but in no event 
less than a reasonable degree of care and warrants that it has in place reasonable safeguards against the unauthorized 
disclosure of Confidential Information. 

3, The Confidential Information provided by the Disclosing Party to the Receiving Party shall not be used 
by Receiving Party to prepare any patent application or to provoke any interference with any patent application that has 
been or may be filed on behalf of the Disclosing Party, including any continuing and/or divisional patent application that 
has been or may be filed. The Confidential Inft>rmation of the Disclosing Party shall not be used by the Receiving Party 
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incombustibility, and is chemically stable (shown in p. 6 
lines 12-19 of the present specification) . Accordingly 
^^ionic liquid" does not refer to a state of matter which 
dependent upon conditions but refers to specific compounds. 
Consequently "Honic liquid" does not mean the ionic 
compounds in the liquid form. 

Therefore, we believe that the term ^^ionic liquid" in 
claims 1 to 6 is appropriate. 

2. With rsspecb to Claim Rajftationa - 35 USC § 102 

Angell (USP5962169) discloses the rubbery 
materials (column 5 line 22-23 etc.).. That is to say, ionic 
liquid is rubbery itself in Angell (column 3 lines 16-18) . 
On the other hand, the present invention is the rubber with 
ionic liquid. Therefore, the present invention is 
different from Angell, 

3. With raapsat to Claim Rejections - 35 USC § 103 

Kitano (USP6810225) discloses ionic compounds, but does 
not disclose the compounds correspondihg to the ionic 
liquid. On the other, hand, the present invention uses the 
''ionic liquid". 

Takashima (USP6458883) also, does not disclose the 
compounds corresponding to the ionic liquid, 

Angell also does not disclose the compounds 
corresponding to the ionic liquid. 

Michot (USP6841304) disclose molten salts and the 
molten salts are used in electrolytic solutes (Title) , but 
does not disclose that the molten salts are added to. a 
rubber. Accordingly, the technical field of Michot is 
different from the technical field of the present invention 
related to rubber composition and rubber member. 

Therefore, the present invention is not obvious from 
Kitano, Takashima, Michot and Angell. 
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Novel Polymer Electrolytes Prepared by 
Copolymerization of Ionic Liquid 
Monomers 

Masahiro Yoshizawa,, Wataru Ogihara and Hiroyuki Ohno'*' 

Department of Biotechnology, Tokyo University of Agriculture and Technology, Koganei, Tokyo 184-8588, Japan 



ABSTRACT 

Ionic liquid monomer couples losre prepared by the 
neutralization of ^l-vinyli7nid&Z.QLe wtth z^iji^/swybnic 
add or 3 -sulfopTapyi aery late. ^These Tomc^ Uqiiid 
'^mmomer co'tiphs taeretisdmC^^^^ at room tempera- 
ture and showed low glass transition terrtperature (Tg) at 
-S3°C and -73''C, respectively. These monomer 
couples were copolymerizea to prepare ion conductive 
polymer nmtrix. Thus prepared ionic licfuxd copolymers 
had no carrier ions, and they showed very low ionic 
conductivity ofhelozu 10~^ S cm"^ Equimolar amount of 
Uibi^ {liTFSI) to 

imidazolium^sWTuhit was men added to generate earner 
ions in the ionic liquid copolymers, Poly(vmylimidazo- 
lium^-co'vinylsulfonate) containing equimolar LiTFSI 



showed the ionic conductivity of 4x10 ^ Scm'^ at 
30°C Advanced copolymer, poly(vinylimidazolium-co- 
S-sulfopropyl acrylate) which has flexible spacer behueen 
the anionic charge and polymer main chain, shaioed the 
ionic conductivity of about 10^^ S an'^ at 30^Q which 
IS 100 times higher than that of copolymer without 
spacer Even an excess amount of LiTFSI was added, the 
ionic conductivity of the copolymer kept this conductiv- 
ity, This tendency is completely different from the topical 
polyether systems. Copyright <S> 2002 John Vmey & 
Sons, Ltd. 



*Con"espondence to; H. Ohno, Department of Biotechnology, 
Tokyo umversity of Agriculture and Tiedmology, Koganez/ 
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INTRODUCTION 



Recently, air- and water-stable (ianicjigm^have. 
been prepared by the combination of imidazolium 
cation and specific organic anions [1^ 2], Since these 
Ionic liquids are composed of only ions, they have 
very hi^h carrier ion concentration and accordingly 
show high ionic co nductivity at room temperature 
[3-7y Furthermore, (iomG liquid'sf are widely spread 
as some kinds of reaction media and extraction 
solvent because of very unique properties, such as 
nonvolatility, low viscosity, and iiTurdscible with 
water when fluorine-rich anions were used [8-11], 

Since ionic liqiuds are empirically known to be 
prepared with imidazolium salts^ several imidazo- 
lium cation derivatives are used to form ionic 
liquids [12]. The ionic liqiuds are generally pre- 
pared through the anion exchange process from 
salts having halide aruons. There are some draw- 
backs such as incomplete anion exchange process 
and contamination of by-products. In addition, this 
method is limited for only several anion species. 
Development of new ionic liquids is not so rapid, in 
spite of unlimited combination of cations and 
anions [1-71- 

We then proposed novel synthesis metliod to 
prepare ionic liquid by the neutralization of the 
tertiary amines with acids (13, 14]. This method is 



Copyright © 2002 John Wiley & Sons, Ltd. 



Received 5 Novmb^r 2001 
Accepted 13 Mfirdi 2002 



19. Dec. 20 0 5 1 8:23 



KURIHARA INTERNATIONAL PATENT 



NO. 9465 P. 5 



590 / Yoshizawa al. 



useful for not only the . improvement of the 
synthesis but also the analysis of the relation 
between ionic liquid structure and their properties 
(glass transition temperature, melting pomt, vis- 
cosity^ ionic conductivity, and so on). This neu- 
tralization method is capable to prepeire desired 
ionic liquid from one pot reaction without any by- 
product. We discussed the effect of coT.mter anions 
on tlie properties of ionic liquid by changing acid 
species [14]. Ionic liquids thus prepared by the 
neutralization method are- suitable as model com- 
pounds for dialkylimidazolium salts. Several im- 
portant factors of amines and acids were obtained 
to design excelleiit ionic liquids. 

Since ionic liquids show excellent ionic con- 
ductivity, fchey are also investigated as a component 
of an electrolyte matrix. Especially, the develop- 
ments of polymer gel electrolytes have been carried 
out with polymers containing ionic Uq^iids [15-17]. 
However, ordinary ionic li<^utds are improper as an 
electrolyte, because component ions of ionic liquid 
migrate along with the potential gradient. This is a 
fetal drawback whexi only target ions such as 
liti\ium cations are aimed to be transported. 

We have reported two methods to solve such a 
problem. One is the zwitterionic-type ionic liquid, 
of which both cation ar\d counter anion were 
tetliered [18]. The otlier is the ionic liquid polymer, 
in which cations or counter anions were attached 
on the vinyl polymer [19], Although zwitterionic- 
type ionic liquids have relatively high Tm, lithium 
transference numbers m these zwitterionlc-type 
ionic liquids exceed 0.5. Such a matrix should be 
suitable for target ion transport when the ion 
migration was accelerated. On the other hand, we 
have designed ionic liquid polymers such as 
polycattdn (20, 21], polycation-type polymer brush 
[22, 23], and polyanlon [24, 25]. T nese are obtained 
by the polynrerization of iomc liquid monomers, 
which have vinyl groups on either cations or 
counter anions. In addition, polymerization of ionic 
liauid is effective to not only transport target ions 
selectively but also generate carrier ions. It is 
widely known that cationic and anioruc polymers 
offer an atmosphere with hi^h dielectric constant, 
Which facilitates the dissociation of lithium salts in 
the matrix [26]. However, polymerization of ionic 
liquid components induced considerable decrease 
in the mobility of ions in the matrix due to the 
increase of glass transition temperatm-e. So far, it is 
necessary to increase mobility of ionic liquid 
structure by introducing spacer between polymer 
main chain and salt structure [22, 23]. 

In this paper, we introduced viiiyl groups in 
both cation and anion, and novel ioruc liquid 
polymers were prepsured by the copolymerization 
of mese ionic liquid monomers. 



CH2=:CH 



EXPERIMENTAL 

Materials 

1-Vinylimidazole and 3-sulf opxopyl acrylate potas- 
sium salt were purchased from Aldricn Chemical 



oA3-CH2CH2CH2--SOr HnCI^Nn^ 

^CHs:cH2 

2 

SCHEA/15 1. Sfructure of ionic liquid monomers obtained 
by neutralfzaffon. 



Co. Sodium vinylsulfonate was purchased from 
Kanto Chemical Co. Inc. These were used as 
received. a,a'-Azobis(isobutyronitriIe) (AIBN) was 
recrystallized from methanol before use. 

Lithium bis(trifluoromethanesulfonyl)imide 
(Li TESI) was a gift from Sumitomo 3M corpora- 
tion. 

Synthesis of ionic liquid monomers 

Alkali metal cations of sodium vinylsulfonate and 
3-9ulfopropyl acrvlate potassium salt were substi- 
tuted into proton by means of cation exchange resin 
(Amberlite IR-120B H AG). Aqueous solutions of 
these acids were slowly mixed with equimolar 
amount of l-vinyUmidazole on ice bath. The 
mixture was stirred at room temperature for a 
day. The aqueous solution was then dned witli a 
rotary evaporator; the residual viscous liquid was . 
pored into dehydrated diethyl ether and stirred for 
1 hotir* The viscous liquid was collected and dned. 
This was repeated twice and the obtained salt was 
dried in vacuo at room temperature. The structure " 
of the obtained ionic liquid monomers as shown in 
Scheme 1 was confirmed with ^H-NMR. 

Polymerization of ionic liqtdd monomers 

Radical polymerization of these monomers was 
initiated with AIBN (1 moI%) in ethanol at 60 *C 
under N2 atmosphere. After polymerization, the 
solution was poured mto a large excess acetone. 
Precipitation was washed with dehydrated metha-r 
nol and dried m vacuo at 60 ''C. 

Methods 

The structure of these lonzc liquid monomers was 
confirmed by ^H-NMR spectroscopy (JEOL a-SOO 
NMR spectrometer). 

The ionic conductivity of the obtained com- 
pounds was measured with the complex-impe- 
aance method using an impedance analyzer 
(Schlumberger Solai*tron 1260 impedance/ gain- 
phase analyzer) with a frequency range trom 
10 Hz to 5 MHz. The dynamic ioruc conductivity 
measixrement system was developed in our labora- 
tory- [27]. All measurements were carried out m a 
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TABLE 1 . Effect of polymerization time on tfie ionic conductivity and Tg 



1 

Flo 
PTb 
PTc 
PT 



not defected 



Polymerization Time Ta ("Q 
(h) 



Monomer 

3 

6 
12 
24 



^83 
-75 



(7/ at SO'^C (S cm' V State 



3,5 X 10-^ 
8.5 X 10-^ 

<io-^ 



Uouid 

Rubber-like 

Solid 

Solid 

Solid 



globe box imder dry Nz atmosphere in the 
temperature range 10 to 60 at a cooling rate 

DSC measurement was carried out with a DSO 
62O0 (Seiko Instruments Inc.) in tihe temperature 
range -150 *C to 200''C at a heating rate of 10 ^^C 



RESULTS AND DISCUSSION 

Ionic liquid monomer 1 was obtained as a yellow 
transparent viscous liquid at room temperature. It 
IS kiiown that imidazolium salts with Suorlnated 
aj-iions, like BF^" CF3SO3" and TFSI", form 
excellent ionic liquid [6]. The monomer 1 was 
obtained as liquid with Tg of -83 ^C. Viscosi^ of 1, 
however, exceeds 1000 cP at room temperature, 
whicli is much liiglier than that for other typical 
ionic liquids. In spite of high viscosity, 1 showed 
relatively high iorxic conauctivit)^ of 3.5 x 10~^ 
S cm'* at room temperature reflecting low Tg. 

Ionic liquid copolymer PI (P implies polymer) 
has been prepared by copoTymerLzation of ionic 
liquid monomer 1 , Table 1 summarizes the effect of 
polymerization tSne on the ionic conductivity and 
Tg for 1. When polymerization time was 3 hourS/ 
rubber-Ske polymer Pla with Tg of — ZS'^C was 
obtained. Only when"" the polymerization time 
exceeded 6 hours, polymeiized 1 became solid or 
powder. An obvious Tg for these polymers cannot 
be detected with DSC apparatus used in. our 
experiments. 

Since the copolymers have no carrier ions 
because both cations and amons were fixed on 
polj^mer main chain, they should show very poor 
ionic conductivity. Pla showed relatively high 
ionic conductivit/ of "8,5 x 10""^ Scm~^ at 30 "^C, 
Tliis can be comprehended to be due to low 
molecular weight matrix, like oligomer. On the 
other hand, the ionic conductivity of Plb, Pic, and 
PI, which were prepared by longer polymerization 
time (> 6 hours), was less than 10"^ S cm^^ at 30°C. 
Although the molecular weight of PI was not 
measured yet, improved average moleculai- weight 
Was confirmed as the decrease of the ionic 
conductivitv. UrJess otherwise stated, polymeriza- 
tion time or PI is 24 hours hereafter. 



In order to chai-acterize PI as an ion conductive 
matnx, three kinds of iTthium salts • (UBF4, 
LiCFaSOa, and LiTFSI) were added equimolarly to 
ixnidazolium salt unit of Pi. Fig. 1 shows tempera- 
ture dependence of the lohic conductivity for PI 
contairung equimolar amount of lithium salt. PI 
witli LiTFSI showed ionic conductivity o? 
7.2 X lO"'^ S cm"^ at 50^C, whereas that of PI with 
LiCFsSOs or LiBF^ was only 10"^ S cm"r This 
considerable mcxease In the ionic conductivity by 
the addition of LiTFSI cannot be explamed only by 
higher degree of dissociation of the added salt. 
TFSn is known to act as a plasticizer [28], and this 
was also cor\£Lrined in the case of this copolymer. 
However DSC measurements were performed for 
these samples, their Tg was not detected clearly. 
The lowenng of Tg cannot be conjQrmed witli the 
DSC measurement, but the increase of the ionic 



-6 



I 

'5 



-7 
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+L1TFSI 





+UCFsS03 



3.2 



3.3 3.4 
1000 rVK"^ 



3.5 



FIGURE 1. Temperature dependence of the ionic 
conductivity for PI containing LiX (X = BF^" CFaSOa", 
TFSI-). [LiX]/[lm+] = 1 .0 (bymol ratio) 
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TABLE 2. Tha ionic conductivity at room temperature and Tg of each monomers and their copolymers 



Tg m 



cri (S cm ') 



Monomer 



Polymer 



Monomer 



Polymer 



1 



-83 
-73 



-31 



3.5 x.]0-^ 
6.5 X 10-^ 



<io-' 

<10-9 



— ; not detected 

Polymers were obtained by the copolymerization of corresponding monohfiers for 24 hours. 



conductivity cleai'ly indicated tt. Similai- tendency 
was reported on the zwitterionic-type ionic liquids 
[18]. 

However PI, mixed with equimolai^ amouiit of 
liTFSI, showed the highest ionic conductivity 
aunong three, this needs further improvement on 
the ionic conductivity. For this, we designed new 
ionic liqmcl copolymer, wWch have flexible spacer 
between vinyl group and counter anion structure 
expecting higher segmental motion. In our pre- 
vious study^ comb-shaped poiyroers were ana- 
lyzed, and iii.troduction of spacer between polymer 
main chain and cation charge was revealed to be 
effective to maintain high ioruc conductivity even 
after polymerization [22. 23],. Table 2 shows Tg and 
the ionic conductivity of ionic liquid monomers 1, 
2, and corresponding their copolymers. Compoun'3 
2 is liquid at room temperature as well as- 1 and 



shows the Tg at -73 ''C. In spite of low Tg;. the lonic 
conductivity of 2 was 6.5 x 10"'^ S cm"^ at room 
temperatttrey which is one order lower than that of 
1. We considered that larger anion structure of ionic 
liquid monomer 2 leads to lower ioruc conductivity. 
Alter polymerization for 24 hours, P2 showed very 
low ionic conductivity (less tlian""lO^^ Scm""^) 
similarly to that of PI PI showed no Tg at the 
temperature range for"this"'meastU'ement, whereas 
P2 showed relatively low Tg of -31 ^C. From these, 
P2 containing littdum salt was expected to show 
higher ioiuc conductivity tlian that of PI. 

Eqtiimolar amount of LiTFSI was also added to 
P2 to compare the ioruc conductivity. P2 containing 
equimolar LiTFSI showed the ionic concluctivity of 
1.2 X 10"^ Scm"^ at 30*'Q which was 30 times 
higher thaxT. that of PI as shown in Fig. 2. This 
diSerence is attributed to the f ollowmg two factors; 




3,3 3.4 3.6 

TOGO rVK^"" 

FIGURE 2. Temperature dependence of \he ionic 
conductivity for PI and P2 containing LiTFSI. [LiTFSI]/ 
[Im*] =* 1 .0 (by mol ratiof 
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FIGURE 3, Effect of LiTFSI concentration on the ionic 
conductivity for PJ^ and P2. 
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one is the increase of free volume bv the introduc- 
tion of spacer §roup, and the other is the improve- 
ment of flexibility of ionic liqxiid moiety by the 
alkyl spacer. The introductioii of spacer was 
effective to maintain the segmental motion liigh. 
The polymer systems having spacers onto cation 
site or both sites aa^e now in progress. 

The effect of LiTFSI concentration on the ionic 
conductivity of PI or P2 was investigated. The ionic 
conductivity increasecfwith increasing the amotmt 
of LiTFSI for both PI and P2 as shown in Pig. 3. The 
ionic liquid copolymers contairung excess LilTSI 
(200 ^ 300 moI%) showed a constant ionic 
conductivity. The ionic conductivity of P2 is 10 ^ 
20 times higher than that of Pi. As mentioned 
above, this difference was based"on the introduc- 
tion of alkyl spacer. P2 malntams relatively high 
ionic conductivity at even Mgh salt concentration. It 
is well known that the ionic conductivity of 
polyether system.s was lowered by the addition of 
only small excess of salts. In poly(ethylene oxxde) 
systems, the ionic conductivity decreased by the 
addition of small excess salts due to the elevation of 
Tg. We therefore attempted to discuss the different 
salt concentration dependence by means of Tg. 
however their Tg could not be deteiinined wi& 
DSC measurement. This problem maybe based on 
the distribution of molecular weight tor these ionic 
liquid polymers. Continuous study has been 
carried out and results will be reported soon. It is 
possible to obtain- more excellent ionic liquid 
polymers with higher ionic conductivit}^ by the 
suitable molecular design of starting ionic- liquid 
monomers. 



CONCLUSION — 

We obtained novel xordc liquid copolymers pre- 
pared by the radical polymerization of 1-vlnylimi- 
dazole neutralized with acids having vinyl group. 
Ionic liquid copolymer showed very low ionic 
conductivity of below 10"^ S cm"^ at SO°C. When 
equimolar amount of LiTFSI was added, .poty- 
(vinylimidazblium-co-vinylsulfonate) showea the 
ionic conductivity of 7.2 x 10"^ S cm~^ at 50 ^C On 
the other hand, ionic licjuid polymer having flexible 
spacer showed the ioruc conductivity of 1.2 x 10"^ 
S cm ^ at 50 ^'C, which was about 20 times higher 
than that for polymer without spacer. It was 
revealed tliat flexible spacer is effective to improve 
the idnic conductivity of ionic liquid polymers after 
addition of salt. 
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Ampiiiphilic ionic liquid dcrivativoa form self-organized 
lamellar liquid crystals wich room temperature ionic liquids. The 
ionic conductivities along the smectic layers have been obtained 
for the samples aligned ia the cells with electrodes. Tho highest 
value is 4,3 X 10""^ S cm"* at 139 'C 



liquid crystals combine order and mobility in molecular and 
supi^olecular levels!*-^ The control of the molecular alignment" 
and the organized structures m tho liquid-crystalliije states can 
lead to the introduction of anisotropic functions for the organic 
optQejecg-oni c and ion-conductive materials. 'J^eccntly, organic 
[ionic iiouidd have attracted attention as electrolytes* and solvent 
media^ for reactions and extractions because of their specific 
propeniea such as ion-conductive, non-volatile, and catalytic 
propertiesj Here, wo intend to obtain amsotropic ion-active 
materials through self-organization of ionip materials. It has 
already been reported that ionic liquid deiivatives having long 
alky] chains alone show thermotropic smectic liquid crystal- 
Unity.^ If these amphiphilic molecules can form the liquid- 
crystalline orgamzcd structures with ionic liquids, their ionic 
assemblies would be expected to funcdon as anisotropic materials 
exhibiting high ionic conductiviiiea. Their physical properties 
may also be easily tuned by due fraction of ionic liquids in the 
mixtures. 

Wc report here a new class of self-organized systems^ which 
is obtained by mixing of imidazolium ionic liquids and its 
amphiphilic derivatives having an octadecyi allcyl chain, as 
shown in Figure L The ionic conductivities of the self-organized 
materials have been e^cammed. 



'CnH2n+i 



JLa;X = BFa, n - 18 2a: X = BF4, n = 2 
lb;X==Pf^e. n = 18 2b:X-PF«, n = 4 

Figure 1. Structure of l-methyl-S-alkylimidazolium salts. 

Two room temperature ionic liquids 2a and 2b, which have 
hydrophilic and hydrophobic counter anions, respectively, have 
been selected to examine the ionic effects on the self-organization 
behavior for the amphiphilic molecules la and lb. These 
compounds are prepared by N-alkylanon of 1-mcthylimidazole 
with 1-bromoalkane, followed by the metathesis reaction using 
HBF4 or HPFfi aqueous solution,^ 

Figure 2 presents a phase diagram for the binary mixture of la 
and 2a as a function of the mole fraction of 2a. Compound la 
alone forms the thermotropic smectic A (Sa) and higher ordered 
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Figure 2. Phase a-ansirion behavior for the mixture of 
la and 2a on the heating runs. 

smectic phases (Sx}« For the mixtures of la and 2a, these smectic 
phases are observed up to the mole fraction of 0.6 for 2a m the 
mixmre. The isouopization temperatui-es are greatly decreased 
with the increase of the fraction of 2a. The mixtures having more 
than the mole fraction of 0.7 of 2a show no liquid crysiallinity , For 
the mixture containing the smaller fraction of la (0.5 wt%), the 
weak translucent gel is obtained below 50 ''C.'' Under polarizmg 
microscopyt opcically anisotropic fibrous aggregadons are 
observed in the isotropic ionic liquid For the mixtures consisting 
of lb and 2b, no smectic liquid^rystalline phase is seen more than 
the mole fraction of 0.4 of 2b. 

lon-conducrive properties have been examined for com- 
pound la alone and the self- organized mixture of la and' 2a, 
which form the smectic liquid-crystalline phases.^ A glass plate 
with comb-shaped gold electrodes (Au) and a pair of Indium tin 
oxide (rrO) electrodes are employed for the amsotropic 
measurements of the xomc conductivity.'^ For the samples 
forming the Sa phase, a homeotropic alignment is achieved on 
the glass surface with gold electrodes, which has been confirmed 
by polailzing microscopy, as schematically illustrated in Figure 
3a. In this case, the iomc conductivities (cTi//) along the direction 
parallel to the smectic-layer plane are measured for the aligned 
samples. In contrast, for Che ITO elecirodeSp the ion-conductive 
values ((Tu.) along only direction perpendicular to the smectic- 
layer plane are not obtained because of polydomain formation, as 
illustrated in Figure 3b. 

Figure 4 shows the ionic conducti vines of compound la and 
the mixture of la and 2a in the molar ratio of 3 ; 2 Cla/2a = 3 /2) 
as a function of temperature. The ionic conductivities for the 
samples homeotropically aligned in the Sa phases are higher than 
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(a) (b) 



Monodomain Polydomaln? 




Figure 3. SchetnaCic illustracions of the smeetic 
liquid-crystalline tonic assemblies aliened wi^ mono- 
domain (a) and polydocains (b) on the glass surface 
with gold electrodes and on the ITO electrodes, 
respectively. Ei Electric field- 
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Figure 4. Ionic conductivities for the smeetic liquid- 
crysialline samples of la alone and la/2a, which ar© 
aligned with monodomain A) and polydomtuns 
(Oi A) on the gloss surface with gold eleca-odes and on 
the ITO elecu'odes, respectively. 

chose for the samples with polydomanis. When the Uquid- 
crystalline order disappears at the isotropization temperatures on 
heating, the ionic conductivities for the samples filled in the both 
cells with Au and ITO electrodes show no difference due to the 
disorder of Che mixture. The formation of anisotropic long-range 
conductive pathways an the homeoiropically aligned Sa phases, 
should contribute to high lomc conductivities. The incoi-poratxon 
of mobile simpler ionic liquid into these pathways enhances the 
ionic conductions in the layered structures- The Ci// values for the 
homeooropically aligned mixture forming the Sa phase become 
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about 10 timris higher than those for la alone, 

It was reponed that the ionic conductivities of isotropic ionic 
liquids alone gradually increase with the increase in 
temperature,^ In contrast, for the mesomorphic materials. in the 
present study, the diseonnnuoua changes of ionic conductivities 
are observed at the phase transidons. For example,, the ionic 
conductivities (cTj^/) for the aligned sample of la increase at 55 
and decrease at 203 ''C on heating. Tliese temperatures corre- 
spond to the Sx-Sa and SA-Isoiropic (Iso) transitions, respec- 
tively. The isotropic ionic conductivity of the single component of 
2a at 30 °C is 2 x IQ-^Scm**',^ while the anisotropic iomc 
conductivity (ai//) of la/2a in die molar ratio of 3 : 2 is 
2 k 10"^ Scm"' at the same temperanue. It is noteworthy that 
the liquid-crystalline assemblies keep relatively high ionic 
conductivities even chough insulating aliphatic layers az« 
incorporated. 

Amphiphilic ionic molecules presented here foim the self- 
organized lamellar liquid crystals with room temperature ionic 
liquids. These liquid crystals are novel funcuonai materials, 
which show anisotropic ion-active properties. The lontc con- 
ductivities of the organized matenais are tunable by the change of 
the fraction of ionic liquids. 
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Chapter 23 

Characterization of Hydrophilic and Hydrophobic 
Ionic Liquids: Alternatives to Volatile Organic 
Compounds for Liquid-Liquid Separations 

Aim E, Visser, W, Matthew Reichert, Richard P. Swatloski, 
Heather B. WiUauer, Jonathan G. Huddleston, and 
Robin D* Rogers'^ 

Department of Chemistry and Center for Green Manufacturing, The 
University of Alabama, Tuscaloosa, AL 35487 



Herein we describe the physical and chemical prop^es for 
several Ionic Liqmds (ILs) comprised of the l-alIcyi-3- 
methylimidazolium cation with the CI", Br% T, BF4', or PFg" 
amons- Liquid-liquid separations with the PF^" ILs are 
highlighted to illustrate dieir use in '"Green Chemistry,'* as 
alternatives to traditional organic solvents in separations. The 
partitioning of organic molecules, combined with XL 
characterization, has facilitated successful metal ion extraction 
with both anionic and molecular extractants- This report 
discusses our cuirent results in the characterization of ILs, 
their use in liquid^liquid extraction, new concepts for metal ion 
separations^ and solid-state analyses. 



Introduction 



) 



The curiosity surrounding Ionic Liquids (ILs) has focused on their unique 
properties, particularly as alternatives to traditional organic solvents in synthesis 
(/) and in liquid-liquid separations from aqueous solutions (2-6). Separation 
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processes employing liquid-liquid extractions are a coxnmon industrial unit 
operation, the characteristics of the extracting phase can be fine tuned by 
selection from the many organic solvents available (7). However, most common 
solvents are volatile organic compounds (VOCs), and their associated chemical 
behavior (in addition to their potential toxicity) classifies many, of them as 
hazardous substances. 



Green Chexnistiy 

Early social perceptions of chemistry and the chemical industry centered on 
ixmovations in medicine, technology, and the overall improvement in the quality 
of life. More recently, with the rapidly expanding chemical industry, those 
perceptions have changed and many perceive the industry as a major contributor 
towards air, land, and water pollution. Each year, the U.S. produces millions of 
tons of pollution and, at the same time, spends tens of billions of dollars 
controlling this pollution in the production and cleanup cycle. Therein is the 
indication that redesigning of chemical production prxx;esses may be more 
effective through targeting source reduction instead of "end of the pipe" waste 
management or pollution control (5>. Ideally, pollution would be controlled or 
minimized at the source, recycling would be done if possible, and any 
unavoidable waste production would be handled in an appropriate manner. 

The chemical industry realizes that their sustained fiiture growth, revenue, 
and scientific development may depend on considering how their production 
practices and technologies impact the health and^afety of their emplpyees and 
the environment Major laws, such as the Clean Air Act and Pollution Prevention 
Act, have been the impetus for many of these changes (5). 

Subsequentty, the idea of "Green Chemistry** was established to promote the 
research, development, and implementation of innovative chemical technologies 
to achieve pollution prevention in a technologically and economically sound 
manner (S). All this would be achieved, from the green chemistry standpoint, 
through technologies that reduce the consumption or production of hazardous 
chemicals during the manufacturing or processing of chemical products. A 
decrease in the generation of hazardous chemicals, along with a diminished 
potential for health and environmental problems, would be the direct, result of a 
new paradigm for chemical manufacturing or processing. 

New Solyent Technologies 

Many chemicals that are routinely used are a necessary part of the 
manufacturing process and yet dangerous to health and the environment. The 
penchant for traditional solvent extraction systems has come to a crossroads as 
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the emphasis for sustainable technology and green chemistry considers the 
overaU environmental impact of both the process and waste streams generated as 
a result. 

We have been investigating DLs and polymw-based Aqupous Biphasic 
Systexxis (ABS) as alternatives to traditional organic solvents in liquid-liquid 
separations. Although both areas have their own specialized literature, several 
reports highlight separations based on ILs {3,20-12) or ABS (13-15) where the 
merits of both systems show their potential for implementadon into industrial 
separations systems. Most notable from a "Green** perspective is the fact that 
both ABS and ILs can be used in place of VOCs in liquid-liquid separations. 
Other work on novel solvent media has shown how supercritical water (16,17), 
supercritical CQ2 (9), solventless phases (IS), and fluorous phases (19) can be 
used in an effon to broaden the scope of possibilities available for more 
environmentally responsible processes. 

Room Temperature Ionic Liquids (RTILs) 

RTILs (or, more generally, ILs), as discussed in this book, are a class of 
novel compounds composed entn-ely of ions that, together, melt at or below 
ambient temperatures (or below 150 for ILs) and may be thought to resemble 
molten ionic melts such as NaCl.at 800 ^C, While both DLs and molten salts are 
composed of ions, the presence of organic cadons in ILs interrupts the crystal 
packing and lowers the melting point Depending on the composition, the 
resulting IL can be hydrophobic or hydrophilic. An inhwent property of many 
ILs is their miniscule volatility and easily manipulated properties, characteristics 
of a unique class of solvents. Literature reports are awash with cations for use in 
these systems, including ammonium (20-22), pyridinium (23,24)^ pyrrolidinium 
(25), isoquinolinium (26), and imidazolium (2,27), each with the possibility for 
attaching various alkyl groups to the ring or quaternary onium cation. Depending 
on the type of cation investigated and the length of die alkyl chain, ttie resulting 
salts may have a melting point above room temperature and several crystal 
structures have been reported (25), 

^ our research, we have focused on RTILs iacorporating the l«alkyl-3- 
metiiylimidazolium ([Q,mim]^ cation, as shown in Figure 1 (2-4). The alkyl 
group is usually an n-alkane and increasing the length of the alkane chain jaffects 
the resulting properties (e.g., viscosity » hydrophobicity, and melting point 
(2,27)). For example, [CsmimllPFg] is a liquid at room temperature with a glass 
transition temperature at -75 ""C (2) and [Cioniini]IPF^] melts at 38 °C (J). RTILs 
composed of several different cations or anions (i.e., multi-co mponent mixtures) 
can be envisaged to further expand the family of Ronic liquidsy 

IL cations and anions can be used to control the water miscibility of the 
resulting CLs and the choice of anion has been used to the greatest effect on 
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ccmtrolling hydrophobicity, l-alkyl-3-iriethylimidazolium salts of PF/ are water 
immis cible; BF4" salts are water miscible depending oti alkyl chain length, and 
tetrahaloaluroinate salts are moisture sensitive, Fluorinated anions such as 
N(S02CF3)2' impart even greater hydrophobicity despite increased expense (27). 



R3 




R3 CH3; Rl mif-aOume 



Figure I, Generic cations for l-^alkyl'S-methylimidazolium RTILs. 



The ionic liquids we have used in our liquid-liquid studies are all liquids 
below 40 although 150 ''C is the adopted upper tenxperature lirait for Monic 
liquid' classification^ RTE- is an arbitrary subdivision of ILs that are liquid at 
room temperature and, thus, facilitate their use at ambient temperatures. Here, 
we will provide an overview of organic partitioning and metal ion extraction in 
RTIL-aqueous systems with particular, focus on the effects of anion sele<;tivity, 
increasing alkyl chain length, and cation substitution. 



Experimental 

The experimental details for the data azid procedures reported here are 
provided in references to physical properties (2), organic solute partitioning (29), 
and metal ion extraction (4-6). Distribution ratios as reported here are defined as 
the ratio of the solute concentration in the RTIL lower phase divided by the 
solute concentration in the aqueous upper phase. 

Thymol blue (i) partitiomng was measuried via UV-Vis spectroscopy while 
the remaining organic solute and metal ion studies reported here utilized radios- 
labeled chemicals. The methods for chemical and physical characterization of 
the ILs are described elsewhere (2). 



k 
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Partitioiiing of Organic Solutes 



Our initial results (29) for aromatic solute partitioning between water and 
[C4znim][PF^] indicated that neutral, hydrophobic, aromatic solutes have an 
affinity for the ionic liquid phase, as shown in Figure 2. For the simple benzene 
derivatives investigated, partitioning to the RTIL phase generally increased in 
accordance with the solutes' 1-octanol-water log F values, a common measure of 
solute hydrophobicity. This deixwnstrated that partitioning could be achieved 
and, in many cases, partitioning can be predicted (and/or modeled) using 
traditional solvent extraction parameter. In this manner, ILs behaved as 
conventional solvents^ displaying a noticeable affinity for aromatic solutes and 
low aliphatic solubility. Once in the RTIL pha^e, volatile organic solutes can be 
removed by distillation or p^aporation, as the RTILs themselves may have 
negligible vapor pressure* Many nonrvolatilcs can be removed through contact 
with supercritical COj (9). The solubility of [CnmimJtPF^] in various cthanol- 
water solutions (30) suggests complex solution behavior with increasing mole 
fractions of traditional solvents. 
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Figure 2. DUtribution ratios ofdrgank solutes for the [C^ndmJlFF^J-water 
hiphasic system; data from (29)^ 

With charged or ionizable solutes, a change m the aqueous phase pH 
resulted in certain ionizable solutes exhibiting pH-dependent partitioning such 
that their affinity for the RTILs decreased upon ionization (29). As shown in 
Figure 3, solute ionization affects the partitioning of aniline and benzoic acid 
and may be sufficient to cause several ord^s of magnitude difference in the 
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paxtitiomng. Thus, for ionizable solutes, a change in aqueous phase pH could be 
used for solute recovery fix>i& the RTIL phase. 
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Figure 3. Distribution ratios for aniline (pKi, = 9^42) and benzoic acid (pK^ = 
4.19} in [C^mimJlPFd-aqueous systems as a Junction of aqueous phase pH. 

Other solutes may have substituent groups that hinder their solubility and 
partitioning in RTIL-water systems. We have studied the partitioning of 
unsubstituted, r-butyl-, and sulfonated-calixarencs in a series of [Cnmiin][PFfi]- 
water systems as a function of aqueous phase pH {32). The hydrophobic nature 
of the RTILs offers an environment suitable for solubilizing the unsubstituted 
caIix-4 and calix-6-arenes, where distribution ratios for both molecules are above 
100 and are relatively unaffected by aqueous phase pH- However, the r-butyl- 
calixarenes were- not soluble in these RTELs. Adding sulfonate groups to the 
calixarenes diminished their affinity for the RTIL phase and distribution ratios 
indicate those molecules prefer the aqueous phase. 

We further investigated the pH dependent partitioning of the cationic dye» 
thymol blue» in these systems (J). As depicted in Figure 4, thymol blue preferred 
the RTIL phase in both the zwitterionic and monoanionic form and partitioned to 
the aqueous phase under basic conditions. Of the RTILs investigated in liquid- 
liquid partitioning experiments, thymol blue distribution ratios increased with 
RTIL hydrophobicity (or Hpophilicity) and were highest in [CgmimlCPFd. Solid- 
liquid separations were performed with [Cioniim][PFfi] to illustrate how 
successive heating, mixing, and cooling, combined with appropriate pH changefs, 
could remove thymol blue from an aqueous phase with no retendon of the dye 
upon recrystallization of [CjomimlCPFfi] (J), 

The behavior of organic solutes in RTIL systems illustrates both the unique 
and unexpected characteristics of RTIL-based separations- With respect to 
organic solutes, traditional descriptors provide a good indication of solute 
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affinity for the RTIL phase. For ionizable molecules, ionizatioii has an apparent 
effect on phase affinity and may provide a facilitated means for stripping. Wheti 
selecting or modifying a RTIL, the resulting hydrophobicity has a significant 
intact on partitioning results. 
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Figure 4. Behavior of thymol blue in [Cnmim][PF^] -aqueous biphasic systems 
as a Junction of aqueous phase pH(3)^ 

Metal loD Extraction 

In liquid-liquid separation systems, the hydratcd nature of most naetal ions 
lowers their affinity for the extracting phase. This is the case in RTIL systems 
where hydrated metal ions do not partition to ILs from water (5,52). Therefore, it 
is necessary to change the hydration environment of the metal ion by either using 
organic ligands (iJ-JJ) which provide a more hydrophobic region around the 
metal or form neutral compounds, or with inorganic anions (id) that form softer 
more extractable anionic conq)iexes with the metal. Ideally^ in a biphasic system, 
the extractant would remain in the hydrophobic phase to ensure the complete 
removal of the metal ions from the aqueous phase* Thus, the challenges in 
adapting new classes of solvents to traditional separations include finding 
extractants which quantitatively partition to the solvent phase and can still 
readily complex target metal ions; or finding conditions under which specific 
metal ion species can be selectively extracted from aqueous streams containing, 
inorganic complexing ions. 



Organic Extractants 

Our initial studies used l-(2-pyridyla20)-napthol (PAN) and l-(2-thiazolyl)- 
2-napthol (TAN) as extractants for metal ions in light of their success in 
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polymer-based ABS {37) and traditional solvent extraction systems (38). These 
are conventional metal extractants widely used in solvent extraction appUcations. 
When the aqueous phase is basic, both molecules are ionized, yet they 
quantitatively partition to [C^mimJCPF^i over the pH range M3. As other reports 
indicate (39), the ionization of the molecules plays an important role in the metal 
ion^xtraction, as shown in Figure 5. The distribution ratios for Fe^, Co^*, and 
Cd show that metal ions could be extractejd from the aqueous phase at basic pH 
and stripped firom the RTIL under acidic conditions. 




4 6 8 10 
Aqueous pH 



Aqueous Phase pH 



Figure 5. Metal ion distribution ratios with 0.1 mMPAN(a) or TAN(b) in 
[C^imJlPFgJ-etgueous systems as a Junction of aqueous phase pH(37). 

Macrocyclic ligands such as crown ethers have been widely used for metal 
ion extraction (,35,40). The structure and cavity size of the crown ether is the 
basis for metal ion selectivity and by attaching alkyl or aromatic ligands to the 
crown, the hydrophobicity of the ligand can be adjusted. Others have used 
dicyclohexano-18-crown-6 as an extractant for Sr^ in [Cnniim][N(S02CF3)J 
RTIL liquid-liquid systems with impressive results (ii). However, the conditions 
for optimal extraction in the RTIL system (i.e., low HNO3 concentrations) were 
unexpected and contrary to what is typically employed in solvent extraction (40). 
We have used 18-crown-6, dicyclohcxano-18-crown-6, and 4,4*-(5')-di'(r- 
butylcyclohexano)-18-crown-6 as extractants for Sr**, Cs*, and Na* in 
[QmiJCPFd-aqueous systems (4). Our results for the partitioning fiom HNO3 
are contrary to traditional solvent extraction behavior in that while the most 
hydrophobic extracting phase usually produces the highest dfetribution ratios, as 
expected, distribution ratios decrease with increasing acid concentrations. 

It was noted, howevar, that in die presence of increasing HNO3 (i.e., 1 - 8 
M), the degradation of PF^" to PO4'" significantly increased the hydrophilicity of 
the RTIL phase and produced a monophasic system. Using A1(N03)3 as the 
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source of the nitrate anion leads to an increase in distribution ratios (Figure 6) 
and precluded the degradation of PFg" observed at high concentrations of HNO3. 
Our results indicate that metal ion partitioning is very complex in RTIL-based 
liquid-liquid systems and other factors such as aqueous phase composition and 
water content of the RTILs have dramatic effects on both the metal ion 
extraction and &e stability of the RULs. 
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Figure 6. Sr^* and Cs* distnbution ratios with OJ M di-^tertbutyUdicyclohexano- 
18-crown'6 in [Cnmim][PF^] -aqueous systems ([C^mim]'' [CtsmimJ* 
[Csmimt ( ^) with increasing aqueous phase concentrations ofAl(N03)s (4). 

Anionic Extractants 

The distribution ratios of certain metal ions may be enhanced in the 
presence of coniplexing anions such as halides orpseudohalides. Metal ions such 
as Hg^* have large formation constants witti halides {41) and their effect has 
been observed on the partitioning of Hg^* in other systenjs (42). We have used 
halides, cyanate, cyanide, and thiocyanate as anionic extractants for metal ion 
partitioning in [C4xnim]|;PF^]-aqueous systems (5) and the results are shown in 
Figure 7, For the halides, the formation constants for the series of Hg-I 
complexes are the highest and decrease to those for Hg-F (41). That, in 
combination with the relative hydrophobicity of the complexes, explains v/hy the 
iodide complexes produce the highest distribution ratios. The results with the 
pseudohalideSt however, suggest a more complex partitidning mechanism since 
the Hg-CN complexes have the highest formation constants (47), yet display the 
lowest distribution ratios, similar to data found in other separations systems (42). 



19. Dec. 2005 18:27 



KURIHARA INTERNATIONAL PATENT 



NO. 9465 P. 21 



298 



9 

E 102 

CO 

e 
9 



§ 
1 



■S 10-2 



10' 
ICP 









; 


u ■ SCN- ^ 


Br- 
OCN- 








' ci- 












• CN- 



0.0 0.1 0.2 0.3 0.4 0.5 0.6 
Aqueous NaX, M 
Figure 7. Hg'* distribution ratios with increasing aqueous halide or 
pseudohalide concentrations in (C.^mmJlPFeJ-'aqueous systems (5). 



Task Specific Ionic Liquids 

Through a coUaboration with the Davis Group at the University of South 
Alabama, Task Specific Ionic Liquids (TSILs) have been utUized for nietal ion 
separations. Here, knoira metal ion Ugating groups are incorporated into the 
cationic .maiety of the TL by tethering to the imidazolium cation (Figure 8). 
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Figure 8. Hg^* and CS* distribution ratios from water to RUL phases 
comprised of 1:1 combinations af[C^ml[PFg] and a TSIL (6). 



We have used, thioether. urea, and thiourea derivatized imidazolium cations 
in conjunction with PFi for the extraction of Hg^* and Cd^* (.6). As shown in 
Figure 8, distribution ratiois are typically higher for Hg^* and a change in the 
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aqueous phase pH has a slight effect on the partitiomng. (Thus, alternate routes 
may be necessary for significant metal ion stripping.) Depending on the 
properties of a specific TSIL, it can be used as the extracting phase or, if the 
TSBL melts above room temperature, in a combination with [C4mim]p^Fd as the 
extracting phase (5). 



Physical Properties 

Utilization of RTILs in solvent extraction technologies requires detailed 
knov'ledge of the solvents' physical properties (e-g., hydrophobicity, viscosity, 
density, surface tension, etc.). Both hydrophobic and hydrophilic RTILs have ' 
their niche, although each type may be uniquely suited for certain applications 
including synthesis, catalysis (43), and separations- We have utilized RTILs 
incorporating either the N(S02CF3)2-, and N(S02CF2CF3)2 ([BETI]) anions 
to produce hydrophobic RTILs suitable for use in place of organic solvents for 
separations (2). (It should be noted that certain combinations (though not all) of 
miidazoUum and pyridinium cations and N(S02GF3)2% N(S02CF2CF3)2", and 
other hydrophobic anions are covered by U. S. Patent 5,827,602 {44).) 

Solvent choice can have a substantial effect on the reaction conditions since 
parameters such as solvent polarity, density, and viscosity impact not only the 
solute environment, but also the associated engineering considerations. With 
traditional organic solvents, associated properties include high solute capacity, 
solute compatibility, and fine-tunability, although it is their volatile nature that 
has raised the most concern and initiated the search for solvent alternatives. 

Table I is a suromary of the properties for water-equilibrat-ed [l-alkyl-3- 
methylimidazoliumjlTFd RTILs, as employed in many of our liquid-liquid 
experiments. Our previous results (2) indicate water content has a profound 
effect on certain properties and, in liquid-liquid separation systems analyses after 
equilibration with water provide the most accurate data needed for engineering 
RTIL process^es under those conditions. The data in Table I is representative for 
hydrophobic RTILs in comparison to die corresponding hydrophilic RTILs listed 
in Tables II and HI. 

Despite the fact that RTILs mcorporating the FPi anion sustain biphasic 
systems, after contact with water they contain a significant, albeit non- 
stoichiometric, amount of water. We have shown that, for example, the water 
content for [C4mim3[PF6] decreases to 590 ppm after a period of drying (2), 
although the resulting [C4mim]|PFd is hygroscopic and rapidly adsorbs 
atmospheric water to return to the saturation point. The RTILs can be 
successively dried by heating imder vacuum. Previous partitioning results show 
that an mcrease in the alkyl chain length produces higher distribution ratios (3) 
for hydrophobic solutes, confimiing the increase in RTIL hydrophobicity and 
lipophilicity. 
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Table L Properties of Water-EquUlbrated CComiinirPF<tl RTILs (2) 

. Parameter 

Water content (ppm) 
Viscosity (cP. 25 ^) 
Density (g/mL) 
Melting point C^C) 
Glass transition temperature C^C) 
Thermal decomposition (^C, onset) 
Surface tension (dyn/cm) 

^ na = transition not observed. 



[C^nim] 


fC^nimJ 


[Cstnim] 






fPF,] 


11700 


8837 


6666 


397 


452 


506 


1.35 


1.24 


1.16 


4 


na' 


na 


-86 


-75 


-75 


360 


390 


374 


49.8 


36.8 


34.2 



[Comim]* (n < 4) combined with halides or BF4" (45) produce hydrophilic 
JLs that have subtly different properties, as shown in Tables II and TTT, 
respectively. In reality, a switch from hydrophiliq to. hydrophobic RTILs is a 
continuum that can be modified by changes in anion type and cation substitution- 
Other cation types, for example PR^^ offer alternatives, especially lipophilic, 
hydrophobic systems which have densities less than 1 which could be used in 
conventional solvent extraction flotation systems. In Tables II and HI, it is 
important to consider that these RTILs were characterized after being *dried' for 
4 - 5 h at 70 under reduced pressure. 



Table Properties of Dried^ HydrophiUc FC«mim][Cl] ILs (2) 

[C^mim] [C^nim] [C^mim] 



Parameter 


raj 


fClJ 




Water content (ppm) 
Viscosity (cP, 25^0 


2200 


1130 


890 




716 


337 


Density (g/naL) 
Melting point (^C) 


1.08. 


1.03 


1.00 


41 


na 


na 


Glass transition temperature (**C) 


na 


-75 


-87 


Thermal decomposition (**C, onset) 


254 


253 


. 243 


Surface tension (dyn/cm) 




42.5 


33.8 



" dried « heating to 70 °C for 4 -5 h while sdniiig, under reduced pressure. 



The properties in Tables I-Zn confirm that compromises may be necessary 
when selecting an optimal IL for a particular application. Although many cation 
and anion combinations are available, it is evident that concessions must be 
made between hydrophobicity, viscosity, expense, etc. The majority of these 
anions are non-coordinating, although the relatively high melting point for 
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[C4inim][Cl] suggescs that the cation-aoion interactions facilitate hydrogen 
bonding. 

Table III. Properties of Dried' HydrophUic [C<tiniiii]CBF4] and 







[C4tnimJ 


Parameter 


[BF^J 




Water content (ppm) 


4530 


1870 


Viscosity (cP,25°C) 


219 


1110 


Density (g/mL) 


1.12 


1.44 


Glass transition temperature (^C) 


-97 


na 


Thermal decomposition (^C, onset) 


403 


265 


Surface tension (dyn/cm) 


46.6 


54.7 


* dried «: heating to 70 for 4 - 5 h while stirring, xmder reduced pressure. 



In general, rheoiogical properties of the RTBLs vary with water content For 
example, increasing the water content usually reduces viscosity. For [Csmixn]^ 
the longer alkyl chain produces a more viscous RTIL while [C4mim3^ .by 
contrast, has the lowest viscosity. Small amounts of chloride remairdng.from the 
synthesis may also act to increase the viscosity (46), 

Density in the ILs described here decreases as the size of the cation 
increases since the mass of CH^ is less than that for an imidazolium ring. In 
addition^ as the cation size increases, surface tension decreases which can be 
attributed to greater charge dispersion over the cation and the associated effect 
on the cation-ai^on interactions- Ordinarily, these properties would be expected 
to change in a similar manner; increasing the cation size would increase both the 
viscosity and surface tension. Thus, viscosity and surface tension, both 
parameters that are highly dependent on intermolecular interactions, appear to be 
determined by different interactions in the E^. 

Thermal properties such as glass transition temperatures and the onset of 
thermal decomposition are essentially independent of the composition for these 
[CnmimlCPF^] RTILs. In general, a class of cations will show glass transition 
temperatures around a similar temperature range (25,47,48). Here, the glass 
transition temperature is around -80 ''C and is close to that for other 
imidazolium-based RTILs (45). A glass transition requires cation and anion 
reorganization to produce a more ordered material. The addition of CH^ groups 
to the alkyl chain may require a minimal amount of addidonal thermal energy to 
produce the transition, hence the similarity in these results- Melting points, on 
the other hand, are not observed for [QmimJITF^] and [C^mimllPFfi], further 
supporting the observed tendency for RTIlLs to supercool. 
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When designing a RTIL, fine-tuning the properties can be achieved by 
changing the cation substituent groups, atuon identity, or by mixing two types of 
RTILs with differing, but defined characteristics. For example, the combination 
of a TSIL with [C4iniwGCPFs] can be osed for Hg** and Cd** extraction (d). 
fiicorporating anions to produce hydrophilic RTILs illustrates a wider range of 
properties, as shown in Table n. The effect of anion identity and, as a result, 
water content on the propwties is obvious by comparison of Tables I-HI. 

Recently, we have prepared a series of RTILs containing 1-aikyl- 
isoquinolirdvim cations, combined with the bis(perfluoroethylsulfonyl)imide 
(N(S02CF2CF3)2', [BETl]) anion, to produce a new class of ILs. We have also 
pursued their characterization and use in separations (2(5). The physical 
properties for [CaisoqJfBETI] are summarized in Table IV. 



Table IV. Properties of Water-Equilibrated rC,isoairBETI7 ILs (2d) 
n Water Content (ppm) Melting Point CO Glass Transition CO 



4 17700 


-62.0 


-85.2 


6 16200 


-77.3 


-84.0 


8 14900 


-68.1 


-79.4 


10 6900 


-59.3 


-77.8 


12 6600 


-51.0 


-75.4 


14 6100 


■49.7 


-66.7 


16. 4700 


-48.6 


-61.8 


18 4100 


-47.2 


-59.3 



As the alkyl chain length ranges from butyl to octadecyl within this series, 
the hydrophobicity shows an order of magnitude increase. [CgisoqJIBETI] and 
[Ci4isoq3[BETI] were used in liquid-liquid partitioning experiments, for a series 
of organic molecules (Figure 9) {26). Increasing, the solute hydrophobicity (as 
denoted by the 1-octanol-water log P value produces increasing distribution 
ratios that are slightly higher in [Ci4isoq][BETI] compared to [C8isoq][BETr|. 
congruent with the water content- In particular, the distribution ratios for 1,2,4- 
trichlorobenzene in [Ci^soqJpBETl] is 1280 (6100 ppm water) compared to 
1130 in [CgisoqJIBETI] (14900 ppm water). Using [C4mim][PP<5] (11700 ppm 
water) as the extracting phase, the distribution ratio is 524 (29). 

Here, the composition of the EL makes significant contributions towards the 
partitioning results. Upon equilibration with water, [Ci4isoq][BETI] contains 
more water than water-equilibrated [C4niim][PFa] (2), although the distribution 
ratios for 1,2,4-trichldrobenzene are much higher than in [C4miin][PF^]. It is not 
yet clear whether the extended aromatic region in [Ci^scq]'^ promotes 
interactions with aromatic solutes or if the observed trends are due to the 
extended hydrophobic environment offered by the longer alkyl chains. 
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Crystal Structures 

Since the scientific commumty has suggested an arbitrary 150 as an 
upper limit on 'Tonic Liquid" melting points, there are many cation and anion 
combinations that produce crystalline materials that can be investigated to 
further understand the interactions in these systems. Parameters such as cation 
substitution and anion identity have been selected to investigate their 
contribution towards the solid state structure of ILs. Initial modeling studies 
focused on basic [l-mcthyl-3-ethylimida2olium][AlCl4] systems and two models 
were proposed for the ionic structures; a stack model in which the anions were 
located between parallel pl^es of cations, or one in which CI" interacts with 
each hydrogen on the ring through hydrogen bonding (49,50). The crystal 
structure for [Cimimiri] shows hydrogen bonding between the anion and the C-2 
ring hydrogen (57). 

For the [Caniim][PF6] ILs, the length of the alkyi chains, and the ring 
substitution affects the solid state structure of the ILs. The crystal structure for 
[C2mim][PF6] has been reported (52) and we have obtained the crystal structure 
for [Ci(vmim][PF6] (52). [CiomimJIPFfi], which melts at 38 and has been used 
in separations at elevated temperatures (i), has a unit cell where the cations are 
oriented in a bUayer arrangement A similar arrangCTient has been observed for 
[Cj2Hiim][PFd (24) where the cation structure also displays a noticeable, and yet 
unexplained, hook in the alkyl chain, 

XLs display essentially *sait-like' crystal packing in three dimensions for 
those cations containing short alkyi chains, while a restricted environment about 
the charged region is present in the bilayer structure for' those with longer alkyi 
chains. The interactions are dominated by ion-ion coulombic interactions in PF6' 



19. Dec. 20(15 18:29 



304 



KURIHARA INTERNATIONAL PATENT 



NO. 9465 P. 27 



m 

m 



ill 
;f# 

t 

I 



i 
11 



if 

(Hi'" 

w 
m 



Structures and change with cation size. Thus, it appears that the key to designing 
iLs lies in understanding weak interactions that arise in addition to the 
coulombic interactions. By necessity, the cations . and anions will orient 
themselves in a salt-like manner, but it is the interactions that result from packing 
frustration that are noticeable in the ILs. 

As shown in Figure 10, the cation structures and lattice arrangements show a 
trend in charged vs* non-charged separations and hydrophobic packing in the 
structure. The crystal structure for the cation and packing diagram for [l-bu1yl- 
2,3-dimethylimidazolium][PFd is shown in 10a where no distinct regions of 
hydrophobicity are obiserved, owing to the length of the alkyi chain* Adding a 
methyl group to the C-2 in place of a hydrogen significantly increases the 
melting point of this IL, considering [C4mim][PFJ exists as a glass under 
practical conditi:ons (2). Figure 10b shows the cation and unit cell for 
[CaomimJCPFfi], As shown in the packing diagram, the bilayer arraiigement 
dominates the structure, leaving a charge-rich region for the cation ring and 
anion to interact Figure 10c shows the cation and lattice arrangement for 
[CamimltPFfi], 




Figure 10. Cation and packing diagrams for three imidazolium-based PF^ ILs. 

The alkyl chain length appears to have a significant influence, on the solid- 
state structures for imidazolium-based ILs and results in two major types of 
orientations; salt-lUce and bilayer. Future solid-state analysis will facilitate the 
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understanding of ILs as to what cation and anion structural paiameters contribute 
to stabili2ing the resulting IL and the nature of the predominant interactions. 

Conclusions 

For separations, challenges lie in understanding the partitiomng results with 
regards to those from more traditional systems. Continued study of organic . 
solute behavior in RTIL-based liquid-liquid separations may facilitate a 
molecular level understanding of the partitioning mechanisms for neutral and 
ionic solutes, ultimately providing a predictive tool for their behavior. Exploring 
the driving , forces for organic solute partitioning will contribute to the 
understanding of metal ion extractants, partitioning mechanisms m RTDL 
systems, and allow incorporation of new concq>ts such as Task Specific Ionic 
Liquids for separations. 

Beyond the applications of ILs in liquid-liquid systems, the assortment of 
available and potential hydrophobic and hydrophiJic RTlLs attests to their easily 
manipulated physical properties. Through careful selection of cations and anions, 
significant or subtle changes can be made to the water content and rheological 
properties of the resulting RTILs. Indeed, crystallographic analyses indicate the 
strucnires of ILs result from several types of interactions which depend on the 
cation substitution and anion type. 

The burgeoning field of IL research is producing exciting results and 
continues to demonstrate the potential for ILs in a variety of applications. 
Especially from the "Green Chemistry" standpoint, the chemical and physical 
properties of ILs, combined with their demonstrated utility as solvent 
replacements, warrants further study and exploration. Thus, tiie number of new 
and exciting exan5>les of IL implementation in separations schemes is limited to 
time ahd the imagination. 
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Conductivity study on ionic liquid/polymer complexes 

Masayoshi Watanabe*,Tomoo Mizumura 

Department cf Ch^misify, Yokohama National University. U6 Takiwadai, ffcdosaya-ku, Yokohama 240, Japan 

Abstract 

Ternary salt matures, consisting of triethylnKsthylammonium benzoaie (TEMAB), lithium acetate (LiOAc) and Hthium 
bis(tnffuort)Tnethylsulfonyl) imide (LiTPSI), form stable molten salts at moderate temperatures at certain compositions A 
typical composition of the molten salts is TEMAB/LiOAc/LiTFSI = 7/2/ 1 in molar ratio, and the glass transition 
temperature, the crystallization temperature and the melting temperature are - 52''C, 30*C and J5^C, respectively Quenched 
salt mixture having the above composition forms a stable supercooled liquid at room temperature, and the ionic conductivity 
IS 10 Scm at 30^C and 10 " Scm'^ at 60*C. Some commercially available polymers such as polyacrylonitrile and 
poly(vtny| butyral) are compatible with the salt mixture to result in film forming polymer electrolytes with 1-2 order(s) of 
magnitude of loss of the molten salt conductivity. 

JCeywords: Polymer electrolyte; Molten salt; Ionic liquid; Ionic conductivity 



1. Introductioii 

A rather new concept of addition of small amounts 
of polymers to real or supercooled ionic liquids has 
been ' reported [1-7] to combine high ionic con- 
ductivity and rubbery compliance in resulting poly- 
mer electrolytes. The high co nductivity and rubbery 
compliance are due to the iotiic liQuids[ and to 
entanglement effects of the polymers, respectively. 
The concept may come from the limitation in . 
conductivity of conventional polyether-based poly- 
mer electrolytes, where, in contrast^ small amounts 
of salts are dissolved in the polymers, and the ionic 
motion is coupled with the segmental motion. The 
limitation is based on the inconsistency of simulta- 
neous increases in both the number of carrier ions 
and their mobility, because the increase in salt 
concentration causes the increase in glass transition 



"Contsponding author. 



temperature (T^) of the host polymers. We have 
explored [1-3] that the addition of small amounts of 
high molecular weight poIy(l-butyM-vinyipyridin- 
ium halide) to room temperature molten salts, con- 
sisting of l-butylpyridinium halide and aluminium 
chloride, yields new polymer electrolytes having 
higher ionic conductivity than 10"^ Scm"* at room 
temperature and rubbery compliance. 

The present work is to extend the new concept to 
lithiura-ion-conducdng polymer electrolytes. AngeU 
et al. have recently reported [4-^7] that so-called 
"polymer-in-salt'* electrolytes, consisting of a super- 
cooled nmixture of lithium salts and a small quantity 
of a polymer like high molecule weight poly- 
(ethylene oxide), give ide^l polymer electrolytes for 
lithium battery application, because of the combina- 
tion of Li ion conduction and rubbeiy compliance. 
Aldiough Tg's of the supercooled mixtures of lithium 
salts lie below ambient temperatures, the mixtures 
are thermodynamically unstable at room temperature, 
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because of melting temperatures {T^) of the mixtures 
. being higher than room temperature [4-7].- Thus, 
crystallization of supercooled mixtures might be, in 
general, a serious drawback of Li*^ conducting 
polymer-in-salt electrolytes in terms of the stability 
of the conductivity [8], The present study is an 
attempt to prepare Li -conducting real ionic liquids 
in the mixtures, consisting of triethyimethylam- 
monium benzoate (TEMAB), lithium acetate 
(LiOAc) and lithium bis(trifluoromethylsul- 
fonyl)imide (LiTFSI). Preliminary results on com- 
patibility of the ionic liquids with some polymers 
and on ionic conductivity of the resulting polymer 
electrolytes are also presented* 

2. Experimental 

LiOAc (Junsei Chemical, mp: 2S0''C) and LiTFSI 
(kindly supplied by M. Gauthier of IKEC, mp: 
were dried under high vacuum at 1 80*^0 for 
30 h before use. TEMAB (Mitsubishi Chemical) was 
recrystallized three times from acetone and dried 
under vacuum at l(fC, mp: 99*'C. Polyacrylonitrile 
(PAN, Aldrich) and poIy(vinyl butryral) (PVBu. 
Sekisui Chemical), . which is acetalized poly(vinyI 
alcohol) by ^-butylaldehyde, were dried under vac- 
uum. 

In an argon filled glove box, weighed amounts of 
LiOAc, LiTFSI and TEMAB were mixed in glass 
vials, heated to a temperature just below decomposi- 
tion temperature of TEMAB, ca, 2G0''C, and 
quenched to room temperature. In the case of 
LiOAc /LiTFSI binary mixtures, the mixtures were 
heated up to a temperature, where homogeneous 
melts were obtained, and quenched to room tempera- 
ture. To a certain composition of the ternary mix- 
tures (TEMAB/LiOAc/LiTFSI== 7/2/1 in molar), 
which formed a room temperature ionic melt, PAN 
and PVBu were complexed in order to obtain film 
forming polymer electrolytes. PAN and PVBu were 
dissolved in N,N-dimethylformamide, and to each 
solution were added the molten salt to yield a 
homogeneous viscous solution, which was cast on a 
glass plate, followed by the complete evaporation of 
the solvent, 

Ionic conductivity of both of the molten salts and 
the polymer electrolytes was determined by complex 



impedance measurements using a computer con- 
trolled HP-4192A LF impedance analyzer. 

For differential scanning calorimetry (DSC) mea- 
surements of the salt naixmres, the quenched samples 
were sealed in Al pans in an argon filled glove box, 
and DSC thermograms were measured by usmg a 
Seiko Instruments DSC 220C at a heating rate of 
lO^'C/min under nitrogen atmosphere. 



3. Result and discussion 

In order to realize real ionic liquids at ambient 
leniperatures (room temperature molten salts), the 
phase diagram for the mixtures of LiOAc and 
LiTFSI was first explored. Fig, 1 shows the phase 
diagram for the LiOAc and LiTFSI mixtures, LiOAc 
is a glass fomaing salt, and LiTFSI is a newly 
developed and high dissociative salt [9]. With the 
addition of LiTFSI to LiOAc, and of the 
mixtures continuously decrease up to ca* 50 mol% of 
LiTFSI. This behaviour is similar to that observed in 
the lithium salt mixtures where a certain lithium salt 
other than LiTFSI is added to LiOAc [4], At 50 
mol% of LiTFSI, and crystallization temperature 
(rj are close to 150^C, whereas is 33°C. Thus, 
the quenched mixture behaves as a supercooled 
liquid at around room temperature. However, vis- 
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Fig. L Phase diagram for LiOAc and LiTFSI binary mixtures. 
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cosity of Che mixture is quite high, and the ionic 
conductivity is only 10^* Scm"' even at lOO^'C. This 
combination was not suitable to be complexed with 
polytners to form polymer electrolytes. 

The second strategy which we adopted to obtain 
stable ionic .liquids is to add TEMAB to the LiOAc/ 
UTFSI mixtures. TEMAB has a relatively low 
(99*'C) in organic ammonium salts and aliphatic 
ammonium salts like" TEMAB is eleccrochemically 
stable up to reduction potenUal of Li"*' ion. Thus, 
TEMAB was used to reduce of the LiOAc/ 
LiTFSI mixtures, serving as a stable ionic solvent for 
the mixtures. Fig. 2 , shows the states of the ternary 
salt mixtures plotted against the compositions. Typi- 
cal DSC thermograms of the salt mixtures are shown 
in Fig. 3, where the numerals in the thermograms 
coxrespond to the composition plots having the same 
numerals in Fig. 2. In Fig. 2, solid circles represent 
the ternary salt compositions that are inhomogeneous 
even at 200"*^ whereas open symbols (circles, 
hexagons and squares) represent the compositions 
that form homogeneous mixtures at 200*'C. The 
compositions represented by the open symbols can 
be divided into three groups, depending on the states 




TEMAB / mol% 



Fig. 2. SiQCcs of ternary mixtures, coHsStsting of TEMAB; LiOAc 
and LiTFSI, plotted against compositions: inhomogeneous at 
20O'*C; O, homogeneous ai 200^C and glassy or highly viscous at 
room temperature (rt> by quenching from 20XfQ\ 0, homogeneous 
a: 200*C and inhomogeneous at n by quenching; homogeneous 
at 200*C and ionic liquids at n by quenching. The compositions 
where numerals are cited in the symbols correspond to the 
thermograms of Fig. 3 having the same numerals. 



355 




temperature / "0 

Fig. 3. DSC thermograms of typical ternary mixmres, TEMAB/ 
LiOAc/LiTPSI: 1, 0.1/0.5/0.4; 2, 0.6/0.1/0.3: 3, 0.7/0.2/0.!; 4, 
0.8/OJ/O.l. 

at room teniperature after quenching from 200^^0. 
Open circles correspond to the compositions that 
form either glasses or highly viscous supercooled 
liquids at room temperature, Open hexagons corre- 
spond to the compositions that result in inhomoge^ 
neous mixtures after quenching. Terbai7 salt mix- 
tures having the compositions indicated by open 
squares form ionic melts at room temperature by 
quenching. These differences are clearly reflected in 
the DSC thermograms shown in Fig. 3. The diagram 
shown in Fig, 2 reveals the following facts. The 
binary compatibilities between TEMAB and LiOAc 
and LiOAc and LiTFSI are rather good, whereas the 
compatibility between TEMAB and LiTFSI is poor, 
especially when LiTFSI composition is high. Reflect- 
ing these binary compatibilities in the ternary mix- 
cures, homogeneous melts at 20O''C (ppen symbols in 
Fig, 2) are obtained when the composition of UTFSI 
is low or when the composition of TEMAB is low. 
For the salt mixtures having low TEMAB and high 
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LiOAc compositions, the quenched salts are glasses 
or highly viscous supercooled liquids, as seen in a 
typically shown DSC thennogram 1 in Fig, 3* The 
states of the mixtures are very close to those of the 
binary mixtures of LiOAc and LiTFSI, as mentioned 
in Fig. L In contrast, for the salt mixtures having 
high TEMAB and low LiTFSI compositions, the 
quenched salts form ionic melts. Typical composi- 
tions of the ionic melts are TEMAB /LiOAc/LiTFSI- 
7/2/1 and 8/1/L DSC results for these two com- 
positions are shown in the thermograms 3 and 4 of 
Fig. 3. of the mixtures are SS^'C and 52*C, 
respectively, and T^'s are well below room tempera- 
ture* The quenched salts are found to form stable 
supercooled liquids at room temperature for at least a 
few days. 

In Fig. 4 h showti temperature dependence of 
ionic conductivity for the salt mixture of TEMAB/ 
LiOAc/UTFSI = 7/2/1, abbreviated as (721) in the 
figure. The temperature dependence was measured 
with descending temperature. The salt mixture was 
kept at each temperature for ca. 1 h for reaching 
thermal equilibrium. The ionic conductivity is 10"'* 
Scm"' at 30*'C and 10"' Scm"' at 60°C. The 
continuously changing conductivity indicates that 
crystallization of the melt does not occur down to 



-2 




.9 I . I . t . f , — I — — , — \ — , — t 
2.S 2.3 3.0 3.2 3.4 3.6 3.8 4.0 

•1000/T/K-1 

Fig. 4. Temperature dependence of ionic conductivity for lomc 
melt (721) (TEMAB /LjOAc/LiTFSI - 7/2/1) and its polymer 
complexes with polymery lonitrile (PAN) and poly(vinyl butyraJ) 
(PVBu). The composition of the polymers are respresented in unit 
mol % in the figure. 



- lO^'C, within die time scale of measurements. The 
salt mixture having the composition of TEMAB/ 
LiOAc/UTFSI= 8/1/1 gave almost the identical 
conductivity^ with the present mixture. 

Compatibility of the ionic melts with polymers has 
been explored by using several commercially avail- 
able polymers. Up to now, PAN and PVBu give 
reasonable compatibility with the ionic melts. Ionic 
conductivity of the polymer electrolytes is shown in 
Fig, 4, The addition of die polymers to the ionic 
melts imparts film forxwing propeny to the resulting 
polymer electrolyte filnis, with a loss of 1-2 order(s) 
of magnitude of the molten slat conductivity. How- 
ever, with standing the polymer electrolytes at room 
temperature for a long time, crystallization and phase 
separation of the salt mixtures were observed. To 
explore other properties than ionic conductivity, such 
as transport number and electrochemically stable 
potential-window, is currently being conducted m 
our group. 
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Molecular Dynamics Study of the Ionic Liquid l'/i-ButyI-3-methylimldazolium 
Hexafluorophosphate 
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We report the results of a molecular dynamics study of the ionic liquid l-rt-buiyJ-3-mRtbyiimidfisoHv^ 
^®5^^?.?PPliosp^ate [bmim][PFd, a widely studied ionic liquid. An aU-aiom force field is devdoped using 
a combination of density ftmctional theory calculations and CHARMM 22 parameter values. Molecular 
dynamics simulations are caxried out in the isothermal-isobaric ensemble at three different temperatures. 
Quandties computed include infrared frequ^ncieSp molar volumes, volume expansivities, isothemial com- 
pressibililiies, aelf-difftisivities, cation-amon exchange rates, rotational dynamics, and radial distribution 
functions. Computed thermodynamic propertiea are in good agreement with available expenmeiital values. 



1. Introduction 



The term l ionic liouldl aL) has been corned in recent years to 
desciibe a class of organic salts that are liquid in their pure 
state at or near room temperature. ^ Some of the more widely 
studied ILs consist of a heterocyclic cation based on a substituted 
pyridine or imidazole and an inorgamc anion. Early studies 
focused on compounds with the [AlCU]" amon, but these liquids 
proved to be unstable in air. This problem was- overcome^ 
through the use of alternative anions such as [BF4]-, [NOa]-, 
[CH3COO]", and [PFe]". Those stable JLs have attracted a great 
deal of interest because of their potentiai as nonvolatile (and 
hence potentially environmencaliy bemgn) solvents.^ They also 
show potential for a range of other applications including 
separations, industrial cleanings fuel ceils, lubricants, and heat 
transfer fluids,'^ 

Compared to conventional organic solvents, relatively little 
18 known about the thermodynamic and transport properties of 
ILs and how these properties relate to chemical constitution and 
structure. In addition to experimental efforts directed at this 
problem, a few theoretical studies have been reported in which 
quantum chemical^ and classical condensed phase simulationfi^' 
were used to examme stnictttral and thennophysical properties 
of water-stable ILs, Other theoretical studios have focused on 
chloroaluminate-based ELs.^*"^^ 

The presecit work reports results of single molecule quantum 
chemical calculations and condensed phase classical molecular 
dynamics (MD) simulations on the ionic liquid l*«-butyl-3- 
methylimidazolium hexafluorophosphate, which will be abbre- 
viated as [bnum][PF6]- Figure 1 shows the optimized geometry 
of a single [bmimJCPF^] molecule in the gas phase,- obtained 
from a density functional theory calculation (see below). The 
atom labels shown in this figure will be referred to throughout 
the work. 

2. Force Field 

The extent to which a classical molecular simulation ac- 
curately predicts theraiophysical properties depends on the 
quality of the force field used to model the interactions in the 

* To ^bom cQirespcn^nce shouM be addressed. E-mail: ed@nd.edu. 
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Figure 1. Geometry optSatiacd structure of tbmimjCPPsJ with atom 
labeifl noted, 

fluid In this work, we have used a standard molecular 
mechanics force field, 1* with ftmctional form 

bonds Anglos 
dihcdiBla impropsi 

J>i 1 L\ '^r; / \ / J ^ai 

where Vrox is the total energy of the system, hamionic potentials 
govexn bond length, bond angle, and improper angle motioA 
about nominal values r^, $0, and v^o and dihedral angles are 
modeled using a standard cosine series. The Lennard- Jones 
parameters for unlike atoms, and rmin,^, are obtained using 
the Lorentz'-Berthelot combinmg rule. Couiombic interactions 
arc modeled using fixed partial charges on each atom center. 
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TABLE 1: Partial Atomic Charges and Lennard -Jones 
Parameters Used in This Work 



acorn 



(«) 


(A) 


(icJmol-0 


acorn 


qt 






C4 


""0,141 


3.600 


0.209 


Hd 


0.055 


2.680 


0.117 






3.OOO 


O.2O9 


• Hio 


O.OOl 


2,680 


0.117 


Ni 


0.1 u 


3,700 


0.837 


C9 


. 0.256 


4-020 


0.234 


C2 




3.600 


0.209 ' 


Hii 


-0.029 


2.680 


0.U7 


Ns 


0.133 


3,700 


0,837 


H12 


-0.099 


2.680 


0,117 


yv 


U.Iol 


2.936 


0.033 


CiD 


-0.209 


4.0S0 


0.326 ■ 


TT 






0.033 


Hl3 


0.051 


2.680 


0,100 


"1 


0.177 




W, 17^ 


£ll4 


U.U*RI 






Gi 


-0.157 


4.550 


0.084 


Hl3 


0.075 


2,680 


0.100 


H* 


0425 


2.540 


0.092 


Pi 


1.458 


4,30 


2,448 


H, 


0.073 


2.640. 


0,092' 


Fi 


-0.421 


3.400 


0.377 


Hj 


0.142 


2.640 


0.092 


Fi 


-0,426 


3.400 


0.377 




0.093 


4.550 


0.084 




-0.368 


3.400 


0.377 




0,055 


^680 


0.092 


Pi 


-0.368 


3.400 


0.377 




0.045 


2,680 


0.092 


Fs 


-0.364 


3,400- 


0.377 


Ci 


r0.122 


4.020 


0.234 


F« 


-0.414 


3.40O 


0.377 



With the exception of n)» and .^o all cation intramolecular 
and Lonnard- Jones parameters la eq 1 were taken directly from 
the CHARMM 22 force field^^ ^g^^ without modification. 
For the anion, bond lengch and bond angle parameters were 
derived from ab initio calculations, as described below. 
CHARMM 22 Lennard-Jones parameters were used for phos- 
phorus and fluonne. 

The minimum-energy geometry of the [bmim] cation and the 
[PFd anion was determmed by poiforroing ab initio geometry 
TABLE 2; Bond, Angle, Dihedral, and Improper Force Constants 
force force conscant 
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optimizations on the isolated cation and anion at the B3LYP/ 
6-3 U+G* level of.theoiy using Gaussian 98, Cation and anion 
charges were set at +1 and —1, respecnvely. The resulting 
minimized structure was used to set ro, So, and Geometry 
opiimizaaon was also carried out on the [bmun][PF^] pair at 
the same level of theory. The amon was initially located well 
away (over 10 A) from the caaon and facing the imidazoliom 
C2-carbon side. The mmmium energy structure agrees very well 
with the results of Meng ec al.^ Partial atomic charges wero 
derivQcl from this geometry using the CHELPG^^ method. A 
listing of all force field parameters is given in Tables 1 and 2. 
We note that the pair calculations yielded total charges on the 
cation and anion of +0.904 and —0.904, respectively. In 
addidon, the anion is slightly polanzed (i.e., the charges on the 
fluorine atoms are not symmetric). For most of the simulations, 
we used these partial charges as a simple approximation of the 
induced polarization of the ions that occur in the liquid. This 
turned out not to be essential, given that the anion exiubited no 
preferential onentation relative to the nearest cation, as deter- 
mined by moniionng the fraction of the time a given fluotuie 
atom was closest to the C2 carbon of the cation. We also 
performed simulations on an anion with symmetric charges 
totaling -0.904 and found essentially no difference in the 
computed properties. 

The force constants k(, and ke for the amon were deternoined 
by performing a vibrational analysis on the geometry-optimizBd 
pau: within Gaussian 98. The normal modes corresponding to 



force 



bonds 

C3.4-N,j 

Ci-Nu 

C4-C5 

Ci-Hi 

Q-aj.6 



220.0 
220,0 

400.0 
400.0 
410.0 



bonds 



force consiam 



1.470 
1.483 
1.382 
1.337 
1.361 
1.078 
1.078 
1.089 



C7— H7,5 

C9-C10 
P-P 



200.0 
222.5 
222.5 
260.3 



1.091 
1.096 
1.093 
1.530 
1.534 
1.530 
1.646 



angles 



kfi (IcJ mor » rad-S) 



angles 



^o(deg) 



C5-C7-N1 
Cj^-Nu-Ca 

N1-CS-C4 
Ni-Ca-Nj 
NU-C2-H1 

N3-C4-Ha 



140,0 
130.0 
30.0 
30,0 
130.0 
130.0 
25.0 
25.0 
25.0 



112,6 
125.9 
109.6 
106.B 
107,2 
109.1 
125.5 
130.8 
122.1 



H7.8-C,-C8 

C?,a-Ce.j— Cp.io 

Hi5,i4,is— Cio— H 
F-P-F 



35.5 
35.5 
33.4 
33.4 
58.4 
34J 
34.6 
35,5 
194,1 



109,3 
107.2 
111.5 
109.5 
1U,6 
106.4 
109.7 
107.6 
90.0 



dihedtftl 


i!^CkJmol-i) 


n 


i5(de8) 


dihedral 




n 


d(deg) 




14,0 


2 


180 


Ha-O-Ns-Cj 


3.0 


2 


180 


N1-C3-C4-N1 


14.0 


2 


ISO 


Nu-C5.4-C4,5-H2J 


3.0 


2 


180 




14.0 


2 


180 


Ni.a-Ca-N^.i^Ci.7 


0.0 


2 


180 


H|-C2-Nu-C,.4 


3.0 


2 


180 


Hi"C,-N3.i-C^j 


0.0 


2 


180 


Hj — C4 "C5 ""Hs 


2.0 


2 


180 


Hw-C4J-K3.j--C6,7 


0.0 


2 


180 


Cij-C,,*-Nu-C7.6 


0,0 


1 


0 


C2-Ni;,-C7.6-H 


0.195 


2 


ISO 



dihedral 


kj( OcJ mor^) 


n 


<3 (deg) 


dihedral 


h( (kJ morO 


n 


c5(deg) 




0.0 


3 


0 


Ni-Ct-Cb-Hmo 


0,0 


3 


0 


Ci-Na.i-C-Ca 


0.1 


3 


180 


C7— Ca— C9— Cio 


0-15 


1 


0 


C,-Ni-C7^C5 


0.2 


4 


0 


H7.8"C7'~Q-"H9.lO 


0.195 


3 


0 


Ni-C7-Ca-Cft 


0.0 


3 


0 


H,C-C8-Q,-H,C 


0,195 


3 


0 


H,C-tC9-Cio-H 


0.16 • 


3 


0 








improper 


JS>(kJmol-irad-2) 




Vo(deg) 


improper 


AV.OcJxnol-»rad-*) 






Hi-Ni-Ns-Ca • 


0,50 




0 ■ 




0:5 




0 


Nu-C4^-C3-Q,7 


0.6 




0 









19. Dec. 2005 18:31 KURIHARA 

Dynanncs Study of [bimiii][PF6] 



INTERNATIONAL PATENT 



die fltretchiBg of a P-P bond and the beading of an F-P-F 
bond were identified using Gausavaew, a jgmphical interface for 
the Gaussian program, and the computed vibrauonal frequencies 
weiB used to calculate the corresponding harmonic force 
constants. The geometry and force constants for the cation and 
anion are summagzed in Table 2. In a similar majoner. 
vibrational assignments for the major, fundamental modes of 
the pair were made apd compared against expcirimenral infirared 
spectroscopy measurements, described in the Results 

section. 

,. 3. Simulation Details 

Molecular dynamics simulations of 300 [bmim] cations and 
300 [PFfi] anions were performed with the program NAMD" 
version 2.4 in a cubic cell with standard periodic boundary 
conditions. The simulations were earned out in the isothennal— 
isobaric (NPT) ensemble using a Nos6-Hoover barostat to 
maintain a pressure of 0,98 atm. The temperature was controlled 
via Langevm dynamics with a damping factor of 5 ps"^ Initial 
configurations were generated by randomly mscrting cation— 
anion pairs into the simulation box. To prevent overlap, an 
insertion was rejected if any of the newly inserted atoms were 
within 0,75 A of any existing atoms. The initial configurations 
were relaxed using a conjugant— gradient energy minimization 
scheme. For most of the simulations, the mitial cell volume was 
chosen to match the experimental density of the state point of 
interest All of the C-H bonds were held rigid using the 
SHAKE^^ algorithm. The r-RESPA^^ multiple time-stepping 
algorithm was used with a time step of 2 fs for bonded and van 
der Waals interactions and 4 fs for electrostatic interactions. 
The dispersion interactions were cut off beyond 20.0 A. A 
switching function, initiated at a distance of 18.5 A, was used 
to bring Che dispersion interactions smoothly to zero at the cutoff 
distance. Long-range electrostatic interactions were computed 
using the particle mesh Ewald method.^^-^ 

Prior to conductmg a production run, the total energy and 
cell volimifis were nionitored until steady-state was reached. 
Equilibration, tttoes varied from 700 to 1000 ps, after which 
production runs lasting 4 ns were started. The thermodynamic 
properties of the system (total energy, pressure, temperature, 
kinetic, and potential energy contributions) were saved to disk 
every 100 time steps, and the atomic coordinates were saved to 
disk every 400 time steps* AH classical simulations of [braim]- 
[PF6] wore conducted at atmospheric pressure (0.98 bar) and at 
the tenq)eratures 298. 323. and 343 K. 

The program NAMD was chosen primarily because it is a 
parallel MD program that scales remarkably well with the 
number of processors, Simulations were run on a cluster of eight 
Dell Precision 530 computois running Linux. Each node coniams 
two Intel Xeon processon operating at 1.7 GHz. A 4 ns 
simulation of 300 IL molecules (9600 atoms) took approxi- 
mately 300 h (12.5 days) to complete. 

4. Results and DlscussSob 

Figure 2 shows the oxperimental*^'^'^ and computed IR spectra 
for [bmim][PFd. The relevant peajfcs in the spectra are labeled 
I-V. The computed frequencies of the major vibrational modes 
agree well with experiment* Using Gaussview, peak I was 
visually identified as the out-of-plane bending of the C2'"N|- 
C5 angle in the imidazolium ring, peak U is the F-P-F bending 
motion, peak m is the in-plane bending of H-C-N angles in 
the imidazolium ring, peak IV is the vibrational motion of 
several atoms in the imidazolium ring, and the peaks in region 
V are the stretching of C— H bonds. The computed frequency 
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Figure 2. Comparison of computed and expeiimencal'^ IR spectra for 

[bmim][PF5]. 

TABLE 3: Comparison of Predicted and Experimental IR 
Spectral Data 



assignment 



pred. freq. 



expL»«fteq. 
(cm-i) 



C2-N1-C3 bend 669.9 622,9 

P-'P-F bend 827.5 815.S 

imidazole H-C-C & H-C-N bend 1 179.5 1 167 J 

imidazole ring bends 1 600.9 1574.6 

aliphatic C—H str. * 3049.2 2939.1 

aliphatic G-H str. 3092.9 2966. 1 

imidazole C-H sor. 3294.8 * 3124.2 

imidazole C-H str, 3299.3 3170.5 




aOO 300 310 320 330 

r#rnp«rattjm (K) 



340 



350 



Figure 3, Companson of predicted and expenmental molar volumes 
of [bmim][PF«] at 0.98 bar. Lines arQ linear fits to the data. 

for peak I occurs at 50 cm~^ less than Che experimental value. 
The agreement betwe^ computed and expezimental frequencies 
for peaks n-IV is excellent. Peaks H and m occur at 15 cm*^ 
less. than the experimental values, and peak TV occurs at 26 
cm"* less dian the experimental value. TTie computed frequen- 
cies above 2800 cm"^ (region V) are overestimated, with errors 
larger than 100 cm"^ (see Table 3). The large errors m the 
predicted H stretching frequencies are likely due to anhar- 
monicity effects,^ 

4,1. Molar Volume. Figure 3 shows the experimental^* and 
computed molar volumes as a function of temperature at 
atmospheric pressure. Results are also given in Table 4. At all 
temperatures, the predicted molar volumes are louver than the 
experimental molar volumes by less than 1%. This level of 
agreement is excellent considering that the calculations are 
purely predictive; none of the force field parameters have been 
adjusted to m^tch experimental data. Table 5 gives the 
breakdown of the potential energy into elecuostatic, van der 
Waals, intramolecular, and kinetic energies. The average 
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TABLE 4{ Comparison of Predicted and Experimental 
Molar Volumes, Voiuzne Expansivities, and Isothermal 

Compressibilities at 0.98 Bar 



V (cm^ mol-Q opxlO^CK-i) . >clO« (bar**) 



cemp (K) 


Sim 




Sim 




Sim 


exp^* 


29%,2 
323.2 
343,2 


207.8 
210.6 
212.9 


208.9 
212.1 
214.7 


5.49 
5.42 
5.36 


6.11 
6.02 
5.95 


36.83 
32.86 
39.20 


41,95 
49.35 
N/A 



electrostatic contributioa to the potential energy at 298 K is 
—256.7 IcT/mol, which is twice as large as the van der Waals 
mieracuons and is 46 kJ/moI larger than the value reported for 
[emim][AlCl4].^ The total intramolecular poienual energy, 
calculated by summing the bond» angle, dihedral, and improper 
angle energies, is 89.6 kJ/mol. This suggests that properties are 
most sensitive to the electrostatic portion of the force field. 

4.Z Volume Expansivity (ap). The volume expansivity 
quantifies the extent to which the volume of a fluid changes 
with temperature at constant pressure, and is defined as 



(2) 



The volume expansivity can. be computed by running a series 
of simulations at the same pressure but different temperatures. 
Because the change of volume with temperature is appxxDxi- 
matcly linear, a? can be calculated using eq 2 by fitting a straight 
line to the simulated molar volume data. Table 4 lists the volume 
expansivities computed using this approach compared with the 
expeximeatal results.^^ In all caseSi the simulations predict 
expansivities that are lower than the experimental values by 
about 0.6 K"V It is observed that die predicted expansivitioa 
decrease slightly with increasing temperature, in good agreement 
with the experiments. 

4,3. Isothermal Compressibility {/Ct). The isothermal com- 
pressibility quantifies the extent to which the volume of a fluid 
changes with pressure ai constant temperature, and is defined 



(3) 



In this work, the isothermal con^ressibility was calculated using 
the following fluctuation formula^ 



NPT 



(4) 



Table 4 lists the computed isothermal compressibilities com- 
pared with experimental values.^*^ The computed isothermal 
compressibilities are consistently lower than the experimental 
values by 12-33%. There is still reasonably good agreement 
with experiment, however, especially considering the well- 
known difficulty in computing pressures from an atomistic 
simulation, as well as the inherent inaccuracy involved in 
computmg derivative quantities using a fluctuation formula of 
the type in eq 4. The ovwall good agreemteit between the 
calculated and experimental" PVT properties for this fluid 
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TABLE 6: Self-Dimiaion Coefficients and Cohesive Energy 
Density of [bmimlCPFj veiaus Temperature 



tenip 
(K) 




(m^ 8-0 


c 

(Jcm-= 


298.2 


9.70 ±4.1 


S.82 db 4.2 


761 


323.2 


15.4 :k 5.8 . 


11.1 ±5.6 


738 


343.2 


12.1 d: 7.3 


10.6 ± 7.6 


727 



3.0 



TABLE 5: Comparison of the Contribatlon of Various 
Terms to the System Total Energy 


2.5 


temp elect. 

(10 (kJ/mol) 


LJ 

(kjr/mol) 


kinetic • 
(kJ/mol) 


intra. 
(kJ/mol) 


2.0 


298.2 -256.72 
323.2 -256.21 
343.2 ^-255.07 


-89.22 
-87.19 

-85.51 


100.50 
108.92 
116.78 


89.56 
96.74 
102.4 


1.0 



0.6 



0.0 
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Figure 4* Center of mass radial distributloa functions of [bmim][PP5l 
at 298 K and 0.98 bar. 

mdicatcs that the force field does a reasonable job of describing 
The liquid state of [bmim][PF6]. The force field was therefore 
used to compute other thermodynamic, structural, and dynailiic 
properties for which cxp^imentai data do not yet exist 

4.4* Cohesive Energy Density. The cohesive energy densxty 
of a liquid is defined as^f 



(5) 



where AV^^^ is the change in internal energy of the liquid upon 
isotherxEial vaporization into the ideal gas state and is the 
molar volume of the liquid. To calculate the cohesive energy 
density, the average internal energy and molar volume of the 
liquid were obtained &om the liquid simulation, and the average 
intemal energy of the ideal gas was obtained by simulating a 
single [braim][PF(5] ion pair at the same temperature as the liquid 
but at zero pressure (i.e., no periodic boundary conditions). The 
cohesive energy density of [bmim][PFd versus temperature is 
given in Table 6, At 298 the cohesive energy density is 761 
J cm"*^. In contrast, the cohesive energy densities of the heavy 
hydrocarbons hexadecane and naphthalene are 268- and 4X0 J 
cm'"3, respectively.^ The extremely high cohesive energy 
density of [binim][PF6] stems mamly 6om electrostatic interac- 
tions, and explains why these liquids have such low volatility. 

4.5. Molecular Structure. To obtain a better understanding 
of the structure of this ionic liquid, various radial distribution 
functions (RDFs or g(r)) were computed at each of the 
statepoints. The center-of-mass RDFs for cation-cafion, cation— 
anion, and anion-amon pairs at 298 K and 0.98 bar are shown 
in Figure 4, ThQ RDFs at 323 and 343 K are qualitatively very 
similar and are not shown here, It is observed that the first 
solvation shell for cation— anion pairs forms at about 4.3 A. 
The second and third cation— anion solvation shells occur at 
10.6 and 17.6 A, respectively. Given the long-range Coulombic 
interactions m this system, it is not surprising to see that weak 
ordering persists beyond 23 A. The first peak in the cation- 
cation RDF occurs at 8 A. and the amon-anion RDF shows, 
two peaks at 6.5 and 9 A. This split peak is the result of 
sequential ordering induced by the cauon-anion pairs. The 
RDFs agree remarkably well with those of Shah et al.,* who 
used a simpler, united atom model for this liquid. 



19. Dec. 2005 18:32 KURIHARA INTERNATIONAL PATENT • NO. 9465 P 40 

jjyaaaucs stmy or tbmimjCPFeJ / Pkys. Chem. B, VoL W6, No. 49, 2002 12811 



S3 



0 

,.0 







1"" 










1 1 

P-C 

P-C 

P-C. 


I ; 

1 




















U 













































































B 10 12 14 18 18 20 
r(A) 



Figure 5. Atom— atom radifll dlstribotioft flincnona for [bmmiJ[PF^] 
at 298 K and 0.98 bar. 

Further insight into the liqmd stnicwre can be gained by 
examining the site— site pair RDPs. RDFs for the phosphorus 
^tom of the anion and different imidazoUum carbon atoms of 
the cation are shown m Figure 5. It is observed that the anion 
prefers to associate with the Ci carbon of the imidazolium ring. 
The C2-H1 atom pair constitutes the largest amount of positive 
charge on the cation, so it is not suipiising that the anion 
interacts more strongly with this part of the cation. In addition, 
the locations of the first peaks in- the RDFs of Figure 5 are 
consistent with the ab initio optimized geometry of the cation— 
anion pair shown in Figure 1. 

The RDFs suggest that the mdividual ions in [bmim][PFfi] 
tend to aggregate into well-deiinBd ion clusters. By Integrating 
g{r) out to the location of the first imiiimum, the coordination 
number for the first solvation shell, N, can be calculated via 
the following equation • 



(6) 



where p is tho bulk density and rsh«u is the first nummum'm 
^(r)., Using rjhcu equal to 8.36 A, the calculated coordination 
number for the cation^ardon first solvation shell is 6,8. 
indicating that each ion is surrounded by a cage of nearly seven 
other counterions. Similarly, the cation-cation cooxdinatioxi 
number is 6,2. and the anion-anion coordination number is 6.1. 
Therefore, the total coordination number for the first solvation 
shell of art ion is nearly 13. Note^ however, that the coordination 
numbers computed from oq 6 are sensitive to the choice of r^^. 
For example, setting rsteu equal to 8.0 A reduces the cation— 
anion, cadon-cation, and. anion— anion coordination numbers 
to 6,2, 5.0, and 4.9, respectively. 

4,6. Diffusion. The coefficient of seIf-di££uslon for a fluid 
or solid can be calculated using the Einstein relation^ 



(7) 



where the quantity m braces is the ensemble-averaged mean 
square displacement (MSD) of die molecules and r,- is the vector 
coordinate of the center of mass of ion {. The mean square 
displacements of the canon and anion at 298 K and 0.98 bar 
are shown in Figure 6. The MSDs are linear up to roughly 1 
ns, after which they. show a slightly nonlinear increase* up to 
about 1.8 ns, followed by a second linear regime. A block 
averaging techmque^ was used to obtain the MSDs^ so the 
results are most reliable ai short times. For this reason, the self- 
diffusion coefficients for the cation and anion vvere obtained 
by fitting the slope of the linear region &om 200 to 1000 ps to 
a straight line and applying eq 7. Results are given in Table 6. 
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Flgnre 6. Center of mass mean aquare displacements of [bmim][FF^ 
at 29B K ai^d 0.98 bar. 

At 298 K and 0.98 bar, the predicted self-diffbsion coefficients 
of the cation and amon are 9,70 x 10"^^ and 8.82 x 10"^^ m^ 
5"^ respectively. For companson, die self-difiusion coefficient 
of water^^ at 298 K and 1 atm is 2,3 x lO"' s'^ This result 
is consistent with the fact that [bmimJDPFg] has a much higher 
viscosity (450 cP)3<> tlwi water. If one assumes that the Stolces- 
Einstem relation can be applied to this system, then the self- 
diffusion coefficisnts for the ionic liquid can be estimaied from 
the relation. 



7lL 



(8) 



where Dil — l/2(Z?(bmim] + ^pPej) and 77 Is the viscosity. Using 
eq 8 and ?7HaO = 0,9 cV?^ the estimated self-diffusion coefficient 
for Che iomc liquid is 4,6 x 10"^* m^ s"^ which is fairly close 
to the computed value, especially considering the simplifying 
assumptions behind the Stokes^Einstein relation. The computed 
self-diffusion coefficients fbr [bmim] and [PF6] arc roughly 10 
times smaller than those reported by de Andrade et al.,^ who 
computed self-diffusivides for l-ethyl-a-methylimidazolim 
([emim]) and [AlCL] ions and are also 10 times smaller than 
the experimental result for [enaim],^^ This is not surprising 
because the [exnim] and [AlCU] ions are smaller than the ions 
m this study and because the reported intermolecular potential 
energy' for [emlm][AlCl4] is not as great as it is in this system. 
The computed self-diffusion coefficient for (FFe] is also seven 
times lower than that calculated by Hanke et al.^ firom a 
simulation of dimethylimidazolium ([dmim]) [PFd at 400 K, 
Much of this difference could be due to the fact that the 
simulations of Hanke et al.^ were at higher temperatures than 
the present work and because [dmim] is a smaller cation than 
[bmim]. We also note that the simulations m both ref 6 and 9 
were over a much shorter time scale (100 ps) compared to the 
relatively long simulation times of this work. 

The self-diffusion coefficients for [bmimj and [PFs] could 
not be detennined reliably by fitting the mean square displace- 
ment data at time^ longer than 1 ns. This was due to Inaccuracies 
in the data caused by insufficient sampling. These inaccuracies 
wore apparent in tho relatively large (but not systematic) 
anisotropy in the diagonal terms of the solf-difftisiviiy tensor 
at times greater than I ns. Long-time anisotropy m a homoge- 
neous system indicates that even longer simulations than those 
in this work are needed to accurately determine MSDs above 1 
ns. One must therefore leave open the possibility that the acmal 
self-diffusivity, which by definition is a long-time quantity, may 
differ from diat which is compiuted over Uiese relatively short 
time scales. Nevertheless, the values reported here appear to 
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Bigure 7. Time dependence of to [C«go(f)] for [bmim][t»F6] at 2W K 
and 0.9.8 bar. 

be reasonable estunatoa of the actual self-diffasivity , on the basis 
of the seuiiquantitanve agreement with the results obtained from 
the Stokes -Einsreiii model. The validitx of the self-diffusivity 
can also be tested through use of the following diffusion model, 
Recalliag that each ion is suxrounded by a cage of other ions, 
. we hypothesize that ion diffusion involves a two-step process. 
In the first stop, neighboring ions are displaced enough to disrupt 
a cage and form a diffusion pathway. In the second step, an ion 
leaves its cage, moving to another neighboring cage site. To 
tost this hypothesis; the characteristic time for the breaJdng of 
ion cages was determined through use of a cage corrolation 
function.33 The cage for a given ion is dwnnined by the 
neighbor list for that ion, defined as 

; (9) 
wbsKfir^ is the Heaviside flmction 

aiici rcut is a cutoff radius for the neighbor list. T^e cage 
correlation function is defined as 

where is the number of ions that have left ion j's original 
neighbor list at time f and ^ is the Heaviside function. In 
computing the neighbor lists for all the ions, rcut was set equal 
to the location of the fkst nmrnnum of the cation/anion radial 
distribution ftinction. Figure 7 shows In CcagcCO for the system 
at 298 K, There is a rapid decay of Cc»,c(0 at short times (r < 
0,5 ns), which is attributed to vibrational motion of ions near 
the boundary of a cage. After the period of imtial rapid decay» 
Ccis^it) plateaus and then decays ejtponentiaily after ^bout 0.8-. 
2.0 ns. It is mteresting to note that 0.8 na is about the same 
time in which a shift in the slopes of the MSDs is observed in 
Figure 6. By fimng an exponential function to the region from 
0.8 to 2.0 ns, a ome constant of 1.53 ns is obtained for the 
decay of the ion cages. This time constant is indicative of the 
average time needed for any particular ion to leave another ion*s 
neighbor list 

To investigate the rate at which anions exit the ion cages, 
Ccagft(0 was computed by considering only the anions surround- 
mg a given canon. A plot of In Cgag^CO at 298 K computed in 
this way is qualitatively similar to Figure 7 and is not shown 
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here. Th^ time constant obtained by fitting the mtcOTediate 
region to an exponential' is 3.46 ± 0.04 ns. Likewise, the rate 
at which cations exit the ion cages can be studied by computing 
In CcAao(0 by considering only the cations surrounding a given 
anion. The time constant computed for cation departure at 298 
K is- 3.25 ± 0.03 ns. It is not surprising to see that the time- 
constant computed ftom Figure 7 is about half the time constants 
for cation and anion departures. This is because whenever a 
cation departs from an amon's neighbor list that cation 
simultaneously sees the amon depart from its own neighbor list. 
Thus, two cages decay whenever one ion moves out of another 
ion's cage. 

The cage correlation results suggest that diffusion m ionic 
liquids can be modeled semiquaniitaiively as an activated 
hopping process. Each ion is assumed to reside on a three- 
dimensional lattice, with lamce apacings A given by the distance 
between pealcs in g(r). Ions execute a random walk between 
sites with a charactenstic time r, assumed to be equal to the 
time constant for cage decomposion. The solf-difftision coef" 
ficient for this model is simply^ 

^sdf=g; (12) 

To compute the self-diffusion coefficient for [braim] usmg this 
model, we set A = 8 A and r = 3.25 ns. The resulting ZP,oif for 
[bmim] is 3.28 x 10-»> m^ s^K Using A 8 A and r = 3.49 
ns for [PFfi] gives a seif-diffusivity of 3.08 x lO'^^ m^ s"* 
These results are roughly three times higher th«T ]i the computed 
self-diffusivitiea and a factor of 10 greater than thai estimated 
from the Stokes— Einstein model. It should be expected that this 
simple model over predicts the seif-diffusivity, because it does 
not account for correlated hopmg motion. Jn the real system, it 
is likely that an ion that leaves a cage will frequently return to 
its original location, thus decreasmg die overall displacement. 
The simplified model does not account for this but rather 
assumes that ions successftiUy thermalize in a new cage after 
each hop. Nevertheless, this simple model appears to capture 
much of the dynamics responsible for diffusion in this system. 

4,7, Test of Sampling. Figure 6 shows that at 298 K the 
average ion moves only about 2 A per ns of simulation time. 
This sluggish dynamics raises a SMious concern over the level 
of phase space sampling m the current simulations, as well as 
the suitability of molecular dynamics for these calculations. As 
a smiplc test to determine if the simulations were smck in a 
local potential energy minimum, throe additional 300 molecule 
simulations were conducted for [bmimJCPFis] at 298 K. Each 
simulation was started from a didfferent izutial configuration and 
at different ixuual densities. Figure 8 shows the molar volumes 
of each of these simulations versus time. AH three mdependent 
simulations converged to the same molar volume as that found 
from the original simulation started near the experimental density 
within approximately 200 ps. This uidicaccs that, although phase 
space sampling is still a concern for this system, it appears that 
the present calculations have been run long enough to properly 
sample the equilibrium state of the system. 

4.8. Other Dynamic Properties* 'The rotational dynamics - 
of the anion were mvestigated by computing the cosme of the 
average angle 6 between the vector from the P atom to the Fi 
atom of an anion at time zero and time r. The decorrelation of 
this angle with time gives insight, into the rotational monon of 
the amon. The rotational time constant was obtained by fitting 
an exponential function to long-time decay of (cos 6). The 
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Figure d. Plot of the cell volumes of three izidepeEident simulatuana 
of Cbmiin][PP^ tit 298 K Qnd 0.98 bar. The horizontal dashed Hne is 
the average volume of the oziginal simulation. 

computed time constant for ax)ion rotation at 298 K is ap- 
proximately 28.8 ps. In a similar manner, the rotodonal dynamics 
of the cation were investigated using a vector normal to the 
imldazolium ring. As expected, the rotadonal time constant is 
much longer at approximately 4.3 ns. This huge separation of 
time scales is Significant; although it appears that the transla- 
tional motion of cations and anions is hig^y conelated and slow* 
their rotational motions occur over vastly different time scales* 

5, Conclusions 

Results of a molecular dynamics simulation of D?mnhl[PFd 
are reported. An all-atom force field for the ionic liquid is 
developed usmg a combination of ab initio calculations and 
CHABMM22 parameters. The agreement between tiie e^tpen- 
mental and computed IR spectra is very good, and the vibrational 
motions associated with various peaks m the experunenial 
spectrum are identified. The agreement between experimental 
and computed values of the volume expansivity and isotheimal 
compressibility are good, and the apeomont between molar 
volumes is excellent. The force field was not adjusted to match 
the experimental data. Discrepancies between the simulation and 
expcrmient may be due. to the use of unoptimized potential 
parameters or the neglect of polarxzability. Liquid structure is 
reported in the form of center-of-mass radial distribution 
functions for cation— cation, cation— anion, and anion— anion 
pairs, as well as site-site RDFs, It is observed that the anion 
tends to orient near the Cz carbon of the cation. Self-diffusion 
coefficients for [braim] and [PFg] are computed fcom the slopes 
of the ccntcr-of-mass moan-square displacements of the cation 
and anion, respectively. The reported seif-difElisivities are 2 
orders of magnitude smaller than the seif-diffusivity of water 
at room temperature. The mechamsm for difftiaion of the ions 
is iQvestigaied via the computadon of cage correlation functions. 
A simple random walk diffusion model based on this time 
constant yields a self-diffusivity that is in fair agreement with 
the .calculationsi as does the Stokes -Einstein estimate based 
on scaling with the viscosity and diSusivity of water. Rotational 
dynamics of the cation and anion are investigated via the 
computadon of a comjmon order parameter. The rotational tijue 
constants are indicative of the very slow rotational dynamics 
of the [bmim] cation but the relatively fast rotational motion of 
the [PFfi] anion. 
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